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INTRODUCTION 


This  was  originally  intended  as  a  brief  intro¬ 
duction  to  the  theory  of  waves  scattered  from  a  randomly 
rough  surface.  Now  it  is  no  longer  so  brief ,  but  still 
introductory. 

An  alternate  title  to  this  report  might  be 
"Selections  from  Bass  and  Puks".  The  majority  of  the 
material  contained  herein  may  be  found  in  their  excellent 
monograph  "Wave  Scattering  from  Statistically  Rough 
Surfaces"1.  I  have  rearranged  some  sections,  to  indicate 
their  relationship  and  shortened  others.  The  section  on 
composite  surfaces  contains  results  due  to  Brown2,  which  are 
not  in  Bass  and  Puks,  with  a  derivation  due  to  Dashen3. 

I  have  attempted  to  include  what  I  feel  are  the 
most  important  forms  of  the  average  field  and  average 
intensity  for  both  large  and  small  Rayleigh  parameters.  Then 
I  have  indicated,  in  the  section  on  composite  surfaces,  one 
way  to  combine  these  forms  to  give  a  result  which  I  hope  is 
valid  over  a  wide  range  of  frequencies,  angles,  and  surface 
statistics. 


Shadowing  corrections  are  included  at  the  appro¬ 
priate  places,  but  no  attempt  to  derive  these  results  has 
been  made.  An  excellent  derivation  due  to  Smith4  is 
recommended  to  the  reader.  Some  suggestions  as  to  possible 
extensions  of  the  use  of  shadowing  are  made,  but  their 
validity  awaits  additional  verification. 


It  should  be  noted  that  this  primer  gives  a  single¬ 
bounce  theory,  describing  a  single  encounter  of  a  wave  with  a 
rough  surface.  The  multiple-bounce  theory,  necessary  for 
example  in  shallow  water,  is  significantly  harder.  Attempts 
at  a  theory  have  been  made,  based  on  a  Green's  function  for  a 
waveguide  with  rough  walls5,  but  little  comparison  with  data 
has  occurred.  A  derivation  of  a  loss  per  bounce  using  WKB 
theory  has  proven  useful  in  computing  loss  due  to  the  rough 
walls6 . 

Currently  at  SAI  a  multiple-bounce  theory  from  the 
moving  ocean  surface  due  to  Dashen  and  Spofford  is  under 
development,  with  promising  results  thus  far.  If  the  current 
good  results  continue,  it  is  very  likely  that  this  theory  can 
be  extended  to  a  multiple-bounce  theory  in  shallow  water, 
including  a  rough  bottom. 

Following  is  a  brief  discussion  of  each  section: 
Section  I  briefly  summarizes  what  everyone  should  know  about 
describing  surfaces  probabilistically.  The  reader  should 
scan  this,  to  fix  notation  at  least.  Various  important 
relations  are  introduced  here,  for  later  reference. 

Section  II  tries  to  relate  Section  I  to  reality,  in 
a  brief  discussion  of  some  of  the  rough  surfaces  of  interest 
to  an  acoustician.  A  point  made  there  perhaps  deserves  to  be 
mentioned  here  also.  The  type  of  description  of  a  random 
surface  used  in  Section  I  certainly  does  not  apply  to  all 
surfaces.  It  seems  to  work  well  for  the  ocean  surface  and  it 
may  be  valid  for  the  ocean  basement.  (Note  that  so  little  is 
known  about  the  basement  power  spectra  that  with  a  little 


ingenuity/  one  should  be  able  to  fit  almost  any  basement 
scattering  data.  Talk  about  curve  fitting!)  However,  this 
description  is  probably  not  valid  for  ice  surfaces.  People 
concerned  with  radar  scattering  may  have  to  cope  with  an  even 
wider  variety  of  surfaces,  buildings,  tall  grasses,  trees, 
each  of  which  may  require  a  different  description. 

Section  III  attempts  a  brief  discussion  of  the 
actual  mechanisms  of  scattering.  A  little  tolerance  is 
perhaps  called  for  on  the  reader's  part,  for  the  author  is 
certainly  no  physicist.  Even  so,  the  rather  simple-minded 
concepts  in  Section  III  have  provided  some  insight  to  the 
author,  and  may  prove  helpful  to  the  reader. 

Section  IV  formulates  basic  solutions  of  the  wave 
equation  and  Hemholtz'  equation.  Again,  this  is  necessary  to 
fix  notation,  define  the  appropriate  Green's  functions, 
reflection  coefficients,  etc. 

Section  V  discusses  scattering  from  a  slightly 
rough  surface,  i.e.,  a  small  Rayleigh  parameter.  The  solu¬ 
tion  given  is  the  first  order  term  in  a  series,  where  the 
Rayleigh  parameter  is  the  expansion  parameter.  The  free 
surface,  is  discussed  first,  followed  by  the  rigid  surface, 
and  then  a  general  interface.  In  each  case,  the  mean  field 
<U>  »  0  is  discussed  first,  then  the  average  intensity 
< } U  j 2  >  •  In  discussing  the  intensity,  a  natural  division 
occurs  when  the  scattering  surface  is  "small"  or  "large". 
Here  "small"  refers  to  a  Fresnel  zone  on  the  surface,  so  that 
across  the  scattering  surface  the  incident  wave  has  little  or 
no  phase  change.  In  the  "large"  case,  where  several  Fresnel 
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zones  are  illuminated,  the  curvature  of  the  wave  across  the 
surface  must  be  taken  into  account.  Typically,  in  the 
"small"  case,  a  scattering  kernel  for  a  surface  "patch"  is 
derived,  and  this  kernel  is  simply  integrated  over  the 
illuminated  surface  in  the  "large"  case.  These  distinctions 
do  not  arise  for  the  mean  field,  because,  as  shown,  the  mean 
field  propagates  only  along  the  specular  path,  relative  to 
the  mean  (horizontal)  plane. 

Section  VI  discusses  surfaces  with  a  large  Rayleigh 
parameter,  and  introduces  the  Kirchoff  or  tangent  plane 
approximation.  The  author  has  attempted  to  make  clear  what 
assumptions,  both  physical  and  mathematical,  are  involved  in 
this  approximation.  Perhaps  it  is  worthwhile  emphasizing 
here,  as  do  Bass  and  Puks,  that  this  approximation  is  not  the 
first  term  in  some  perturbation  expansion,  but  simply  an  ad 
hoc  approximation.  Accordingly,  no  error  estimate  is 
possible,  the  only  question  can  be  "Does  it  work?".  The 
solution  does  work  in  many  cases,  of  t  course,  but  these 
remarks  are  worth  keeping  in  mind. 

In  Section  VI  shadowing  corrections  are  introduced, 
as  seems  appropriate  for  very  rough  surfaces.  The  scope  and 
limitations  of  existing  theory  are  indicated. 

Composite  surfaces  are  introduced  in  Section  VII. 
Expressions  for  the  average  intensity  are  derived,  which  are 
simple,  and,  one  hopes,  valid  for  a  wide  range  of  Rayleigh 
parameters.  Shadowing  corrections  are  indicated  as  needed. 


Finally  Section  VIII  considers  the  power  spectrum 
of  a  moving  rough  surface,  specifically  the  ocean  surface, 
rather  briefly.  The  earlier  results  are  extended  to  this 
case  in  a  fairly  straightforward  manner. 
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where  o  is  the  standard  deviation  or  r.m.s  height  of  the 
surface  around  the  mean  value  0,  a2«<C2>. 

This  assumption  is  quite  often  correct,  in  that  the 
sea  surface  has  a  nearly  Gaussian  height  distribution8.  On 
other  occasions,  a  Gaussian  height  distribution  is  assumed 
for  want  of  a  better.  Here  the  central  limit  theorem,  and 
the  mathematical  conveniences  of  normal  distributions  are  two 
motivating  factors  for  this  assumption. 

Two  such  mathematical  conveniences  are  as  follows: 

1)  Linear  combinations  of  normal  random  variables 

are  again  normal.  This  means,  for  example, 

that  a  Gaussian  surface  may  be  written  as  the 
sum  of  two  Gaussian  surfaces. 

2)  Since  differentiation  also  is  a  linear  opera¬ 
tion,  it  follows  that  the  slope  density 

r 

function  for  the  random  variables  and 

&  c 

yy=^y  is  also  normal.  For  a  fixed  surface  in  a 
suitable  coordinate  system, 

w2(vV  "  7*rV  exp\i^Ti  +  "fin  '  (K2) 

*  y  '  lx  y 

2  2 

where  rx2  *  <-57”>r  Fy2  “  <’Sy”>'  * <  >  denotes 

a  ?s  r 

expected  value);  and  *  0  in  these 

^  a  r2 

coordinates.  For  a  moving  surface,  <(^_ )>, 

<TT  T^>  and  <4-r  also  enter  into  the 

dtdx  dtdy 

covariance  matrix  in  the  exponent. 
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Here  the  axes  have  been  chosen  to  eliminate  the 
spatial  cross-term  in  the  quadratic  form  in  the  exponent  of 
W2*  Because  of  the  special  nature  of  the  time  coordinate, 
the  space-time  cross  terms  cannot  be  eliminated. 

If  the  surface  is  isotropic,  r  =r  =r,  and 

x  y 


w2'vV  ■  7^  t-  7^  <V+V}) 


2T4 


(1.3) 


An  additional  important  descriptor  of  a  random 
surface  is  the  correlation  function  W,  given  by 


w(x1 ,y 1 ,t 1 ;x2 ,y2 , t 2 ) 


)c(*2,y2,t2)  >  . 


If  the  surface  is  homogeneous  and  stationary,  then 
W  depends  only  on  the  differences  of  the  arguments  above, 
i.e. , 

w(x,y,t)  *  ~  <  C(x1 ,y1 ,t)c(x1+x,y1+y,t+x)  >.  (1.4) 

Note  that  in  this  definition,  W  is  normalized  so 
that  W(0,0,0)  =  1. 

An  additional  important  density  function  is  the  two 
point  density  function  w(2)  (zj,z2;r,T)  defined  by 

P{zx  <  C(ri,t)  <  zx+dzx,  and 
z2  <  C(rx+r,t+x)  <  z2+dz2) 

2  w*2)  (z1,z2;r,'t)dz1dz2.  (r»(x,y),  ri*(xx,yx)) 
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For  Gaussian  surfaces,  the  following  relation 
between  w(2)  and  the  correlation  function  W  is  valid: 


w 


(2) 


2  no' 


X-W 


2(r,x) 


2W(r,x)z1z2+z22 
2a2  ( 1-W2 (r , x) ) 


]. 


(1.5) 


An  alternative  description  of  the  random  surface 
may  be  given  in  terms  of  the  characteristic  functions,  which 
are  the  Fourier  transforms  of  the  height  probability  density 
functions: 


+  00 


,  ikz  iC(r , x) 

f 1 ( k)  ■  /  w | ( z )  e1K  dz  *  <e  >, 


(1.6) 


+0D 


f2(ki,k2;r,x)  »  //  w(2)(z1,z2;r,x) 


-  ,  i(k1z1+k2z2) 


dz idz2 


(1.7) 


i(k.C(r. ,t)+k9C(r o+rrt+T) ) 

<  e  c  1  >. 


that 


It  follows  from  the  relation  between  w(2)  and  W 


<  C(r1#t)  C(r1+r,t+x)  >  »  a2 W(r,x) 


(1.8) 


&2f ,(k. ,k,;r,x) 


A  further  relation,  to  be  used  subsequently,  is: 


2 

f2(k1  ,k2;r,  x)  -  exp  [^-(k12+k22+2W(r,x)k1k2)  ]  (1.9) 

(an  identity  valid  for  homogeneous  and  isotropic  Gaussian 
surfaces ) . 

The  spectral  density,  or  power  spectrum  F(k  ,k  ,w) 

x  y 

is  given  by 


W(x,y,x) 


—  Re 

«2 


+* 

// 


+<*> 

dk  dk  / 

x  y  _» 


du)F(kx,ky,u>)e 


i  (k»r-u>T) 


(1.10) 


£  -  (kx,ky),  r  =  (x,y). 

For  a  fixed  (homogeneous)  surface,  since  W(r)  = 
W(-r),  it  follows  that  F(  -ic )  *F(  ic )  .  Furthermore  F  may  be 
taken  as  real,  and  by  the  Wiener-Khinchine  theorem,  F>0. 
However,  for  a  moving  surface,  F  is  not  necessarily  even  in 


For  the  moving  ocean  surface,  the  first  order  dis¬ 
persion  relation  for  deep  water  gravity  waves9  implies  that 

F(£,w)  -  F^jo  6(w-u>+),  (1.11) 

where  F^  is  the  first  order  surface  spectrum,  and 
<*>  «  sgn(k  )/gJc  ,  k2«k  2+k  2,  (g-9.81  m/sec2). 


Now  let  Jc*  (  k  ,  k  )  *k  (  cos  4> ,  s  in  4> )  ,  and  r*(x,y) 
x  y 

*  r (cose, sinS ) .  Then  (1.12)  or  (1.14)  becomes 

ao  2  71 

W(r,x)  *  — i-  /  dk  /  kF.  (k,  $)cos(krcos(  6-<|>)-u)  x)d<|>,  (1.16) 

o2  0  0 


and  (1.13)  becomes 

co  2 1C 

W(r)  =  —  /  dk  /  kF.  (k,4>)cos[krcos(  ©— 4>)  )d4>.  (1.17) 

a2  0  O' 

If  the  spectrum  is  separable,  i.e.,  Fj(k,4)« 

F]  (k)H(  ) ,  then  the  (^-integration  may  be  possible.  In 

particular,  if  the  surface  is  isotropic,  then  H*1.  Equations 
1.16  and  1.17  reduce  to 


W(r,x)  =  ±1  J  kF.  (k) J  (kr)  cos(u>,x)dk 
a2  0  1  °  + 


and 


W(?)»  /  kF. (k) J  (kr)dk. 

a2  0  1  ° 


(1.18) 


(1.19) 


If  the  surface  is  cylindrical  or  one¬ 
dimensional,  that  is  independent  of  y,  then  in  (1.11)  F^  ( Sc ) * 
6 (ky ) F 1 (kx) 6(w-w+) ,  where  now  w+*sgn (kx ) /g | kx | .  From  (1.12) 

i  +“ 

W(x,x)  ■  /  F.  (k  )  cos ( xk  -u>  x )dk  ,  (1.20) 

„2  1  x  x  +  x 


and  from  (1.13), 


Just  as  a2=<C2>  is  an  important  parameter  in 
discussing  the  height  distribution,  the  notion  of  a  corre¬ 
lation  time  ?0  and  a  correlation  length  JtQ  are  important 
in  discussing  the  correlation  function. 

Physically  these  describe  a  time  (or  length) 
interval  over  which  different  surface  elements  are  corre¬ 
lated.  On  this  scale  the  interaction  of  an  incident  wave 
with  these  elements  should  be  considered  coherent.  Con¬ 
versely,  over  several  correlation  times  or  distances,  the 
surface  elements  are  uncorrelated,  and  interaction  with  an 
incident  wave  is  expected  to  be  incoherent. 

Naturally  the  transition  from  correlated  to  uncor¬ 
related  is  fuzzy,  as  is  the  definition  of  correlation  times 
or  lengths.  A  local  definition,  emphasizing  small-scale 
roughness,  is  given  by 


V2  W(0) 
r 


(1.22) 


A  global  definition,  involving  all  roughness 
scales,  is  given  by 


1U2  «  ff  W(x,y)dxdy,  JtQ  -  /  W(x)dx. 


(1.23) 


*  One  dimensional  form 


These  varying  definitions  essentially  coincide  if  W 
is  a  Gaussian  correlation  function,  i.e.,  in  one  dimension 


W(x,x)  2  expC—  -  — ). 


(1.24) 


For  more  general  correlation  functions,  they  need  not  agree. 
(Indeed,  the  local  version  may  not  exist,  if  W  does  not  have 
second  derivatives  at  0.) 


For  a  fixed  surface,  the  dispersion  of  the 
derivatives  is  related  to  the  height  dispersion,  and 

correlation  lengths  by 


2  -  02iO! 


t-2(o,0)|,  r  2  -  o,o>|. 

ftv2  j  a*/2 


(1.25) 


For  a  general  surface,  a  knowledge  of  w^  and  W  does  not 
determine  the  higher-order  statistics  of  the  surface. 
However,  for  a  Gaussian  surface,  a  knowledge  of  a2  and  W  does 
completely  specify  all  the  higher  moments  of  the  surface, 
which  is  another  motivation  for  preferring  to  deal  with 
Gaussian  surfaces. 


A  possible  point  of  confusion  should  be  mentioned 
here.  When  speaking  of  a  Gaussian  surface,  what  is  meant  is 
that  the  height  probability  density  function  is  Gaussian  (and 
therefore  so  is  the  slope  density  function).  The  correlation 
function  for  a  Gaussian  surface  may  or  may  not  be  Gaussian, 
(and  the  ocean  surface,  for  example,  does  not  have  a  Gaussian 
correlation  function). 


II.  ACTUAL  ROUGH  SURFACES 

For  the  moving  ocean  surface,  the  Gaussian 
assumption  is  generally  valid.  Further,  there  are  many  semi- 
empirical  determinations  of  the  power  spectrum,  (Pierson, 
Phillips,  Pierson-Moskowitz) 1 1 ,  so  some  degree  of  confidence 
is  possible  when  describing  the  ocean  surface  statistically. 
See  Tables  1  and  2  for  typical  parameter  values. 

When  dealing  with  a  planetary  surface,  the  ocean 
basement,  or  lunar  surface,  the  Gaussian  density  assumption 
need  not  be  valid.  Considerable  disagreement  exists  as  to  an 
appropriate  choice  of  power  spectrum  or  correlation  function. 

Many  of  the  scientists  studying  radar  backscat- 
tering  from  the  lunar  surface  assume  a  correlation  function 
of  the  form12 


W(x,y) 


.  (w(x )  -  e-a  I*  I). 


(2.1) 


However,  it  has  been  shown  that  a  correlation  function  for  a 

3W  3w 

reasonably  smooth  surface  must  have  ■^(0,0)®-^(0,0)»0,  and 

also  have  second  derivatives  at  (0,0) 13.  (The  non-existence 
of  v2w(0)  implies  that  the  surface  has  vertical  faces.) 
Clearly  the  form  for  W  given  above  fails  to  meet  these 
requirements.  Nevertheless,  such  a  function  is  still  being 
used,  primarily  because  with  the  scattering  model  used,  it  is 
often  possible  to  obtain  a  reasonably  good  fit  to  the  data  by 
a  suitable  choice  of  the  parameter  a. 


It  is  quite  possible  that  the  scattering  model 
used  in  lunar  studies  is  inadequate,  which  perhaps  explains 
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why  the  above  correlation  function  is  used  to  fit  the 
data. 


A  kinder  explanation,  and  perhaps  a  generally  true 
statement,  is  that  while  near  0,  the  correlation  function 
does  have  the  required  form  (locally  a  Gaussian  correlation 
function) 


r  2  r  2 

W(x,y)  5  1  -  l  [-5-  x2  +  y2]  ,  (2.2) 

a2  a2 

at  larger  distances,  W  may  decay  like  e~a  ^  rather  than 

e-‘  t  I2. 


Physically  this  implies  that  for  such  a  surface, 
the  correlation  length  is  determined  by  higher  derivatives  of 
the  surface,  curvature,  etc.,  rather  than  depending  only  on 
r.m.s.  height  and  slope,  as  does  a  Gaussian  correlation 
function. 

Using  the  relation  between  the  power  spectrum  and  W 
derived  above,  the  correlation  length  using  the  first 
(derivative)  definition  can  also  be  expressed  in  terms  of  the 
power  spectrum,  e.g.,  for  a  one-dimensional  surface. 


-J 

a*  — od 


pi<W 


dk 


(2.3) 


and  the  theoretical  requirement  that  W''(0)  exist  implies 
that  kx2  Fi(kx)  must  be  integrable. 

An  easy  way  to  meet  this  requirement  in  a  model  is 
to  assume  that  there  is  some  cut-off  wave  number  kc  such 
that 


F1<V  »  0,  |kx|>kc,  (or  F1  (r)  *  0,  |r|>kc).  (2.4) 


Such  an  assumption  implies  that  there  is  some  lower 
limit  to  the  wave  lengths  present  in  a  surface,  or  at  least  a 
lower  limit  to  the  wave  lengths  which  effectively  interact 
with  an  incident  wave. 

This  is  probably  correct  for  the  ocean  surface,  A 
popular  power  spectrum  for  gravity  waves  is  the  Phillips 
spectrum,  which  in  one-dimension,  s  kx"3.  But  there  is  a 
lower  limft  on  the  wave  length  of  gravity  induced  waves, 
below  which  capillary  waves  are  present,  and  the  power 
spectrum  for  capillary  waves  decays  more  rapidly  than  that 
for  gravity  waves.  Capillary  waves  probably  do  not  them¬ 
selves  interact  significantly  with  an  incident  sound  wave, 
but  provide  an  energy  transfer  mechanism  from  the  wind  above 
the  surface  to  the  gravity  waves. 

A  surface  for  which  the  integrability  requirement 
on  kx2Pi(kx)  is  not  met  is  that  of  ice14.  The  power 
spectrum  for  ice  does  apparently  decay  like  kx~3,  and  this 
holds  true  on  very  small  wave  length  scales.  Accordingly,  a 
notion  such  as  correlation  length  has  no  meaning  for  an  ice 
surface.  Such  a  surface,  (within  the  class  of  all  surfaces), 
is  on  the  boundary  of  fractals,  or  diffractals,  and  may  well 
require  somewhat  different  techniques  to  develop  a  scattering 
theory. 

Ice  surfaces,  with  many  small  fractures  and  nearly 
perpendicular  faces,  are  quite  different  from  planetary 
surfaces,  where  weathering  presumably  has  rounded  the  edges, 
and  filled  in  tiny  fracture  zones.  Whether  or  not  such 
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processes,  in  one  form  or  another,  have  occurred  for  the 
lunar  surface,  or  the  ocean  basement,  is  not  clear. 


There  is  no  general  consensus  on  a  "typical"  power 
spectrum  for  the  ocean  basement,  and  some  oceanographers  are 
reluctant  to  concede  that  such  a  notion  has  any  meaning,  due 
to  the  presence  of  sea  mounts,  significant  trenches  and 
ridges,  etc.  Opposed  to  this  is  the  observation  that  data 
(e.g.,  propagation  loss  as  a  function  of  range  and  frequency) 
collected  over  several  areas  of  the  Pacific,  (where  a  thin 
sediment  presumably  allows  interaction  with  the  ocean  base¬ 
ment)  exhibits  a  uniformity  which  suggests  the  same  mech¬ 
anisms  are  operating  everywhere.  If  the  scattering  depends 
on  the  power  spectrum  of  the  ocean  basement,  then  presumably 
there  is  a  typical  power  spectrum. 

An  additional  complication  encountered  when  discus¬ 
sing  bottom  interaction  in  thin  sediment  areas  of  the  ocean, 
such  as  -the  Pacific,  is  the  water  sediment  interface.  In 
regions  where  the  sediment  is  on  the  order  of  30  meters 
thick,  apparently  the  sediment  is  "draped"  over  the  rough 
basement,  and  follows  a  somewhat  smoothed  version  of  the 
basement  contours. 

So  if  it  is  necessary  to  treat  the  basement  as  a 
randomly  rough  surface,  then  presumably  the  sediment  also  is 
randomly  rough.  This  introduces  significant  complications 
into  a  scattering  model,  as  now  a  *K>und  wave  interacting  with 
the  basement  before  returning  to  the  water  column  encounters 
three  rough  surfaces  enroute,  (with  two  transmission- 
reflections,  and  one  reflection). 

As  the  scattering  theory  to  be  described  here  is  a 
single  bounce  theory,  three  encounters  with  rough  surfaces 
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presents  difficulties.  Even  if  the  present  theory  could  be 
extended,  and  for  example,  three  scattering  kernels  obtained, 
whose  iteration  presumably  would  represent  the  scattering 
effects,  the  numerical  difficulties  would  in  turn  become 
formidable. 

One  possible  way  to  avoid  this  difficulty  is  to 
model  the  sediment  as  a  fluid,  with  a  constant  gradient  in 
the  sound  speed,  (and  frequency  dependent  attenuation)  but 
with  an  impedance  match  at  the  water  sediment  interface. 
Then  the  transmission  coefficient  is  1,  and  the  distortion 
produced  by  refraction  at  the  rough  sediment  surface  can 
presumably  be  lumped  into  the  scattering  produced  when  the 
sound  wave  reaches  the  basement. 

The  result  of  this  .assumption  is  a  single  bounce 
theory,  at  least  for  incident  waves  approaching  the  sediment 
above  a  certain  critical  grazing  angle.  For  waves  with  a 
grazing  angle  below  -this,  the  equivalent  ray  turns  above  the 
basement,  and  presumably  the  basement  has  little  or  no  effect 
on  such  waves. 

One  should  also  note  that  if  the  model  allowed  an 
impedance  discontinuity  at  the  water  sediment  interface,  many 
additional  phenomena  would  be  encountered.  Among  these  could 
be  a  combination  of  parameters  producing  a  critical  angle 
below  which  near  total  reflection  occurs.  Such  a  critical 
angle  is  well-defined  for  planar  interfaces.  It  is  con¬ 
siderably  more  difficult  to  introduce  such  notions  when  the 
interface  is  rough,  and  the  notion  of  "grazing  angle"  now 
must  be  interpreted  in  some  average  sense. 


In  thick  sediment  areas,  such  as  the  North 
Atlantic,  with  thicknesses  h  200  meters,  it  is  very  difficult 
for  an  incident  wave  to  reach  the  basement,. as  almost  all  the 
rays  turn  above  the  basement.  And  as  attenuation  in  the 
sediment  is  higher  than  in  the  water  column,  such  waves  are 
essentially  damped  before  they  can  return  to  the  water 
column.  So  the  basement  interaction  is  not  significant. 

Further,  in  thick  sediment  areas,  the  water  sedi¬ 
ment  interface  is  gently  undulating  so  that  either  the  inter¬ 
face  can  be  treated  as  planar,  or  perhaps  as  only  a  slightly 
rough  surface.  Much  of  the  data  from  this  area  can  be 
satisfactorly  modeled  without  recourse  to  rough  surface 
theory15,  so  there  is  little  need  for  a  scattering  model  in 
thick  sediment  areas. 


Of  course,  how  one  determines  that  a  thick  sediment 
model  is  no  longer  adequate,  and  a  thin  sediment  scattering 
model  allowing  for  interaction  with  a  rough  basement  is 
needed,  is  not  very  clear.  One  can  conceive  of  a  sediment  of 
an  intermediate  thickness,  such  that  at  longer  ranges,  the 
sound  waves  refract  in  the  sediment,  and  turn  above  the  base¬ 
ment,  whereas  at  shorter  ranges,  a  significant  fraction  of 
the  sound  energy  in  the  sediment  scatters  off  the  basement. 


III.  DISCUSSION  OF  THE  PHYSICS  OF  SCATTERING 

Before  proceeding  to  a  discussion  of  the  various 
mathematical  models  for  scattering,  it  is  useful  to  examine 
the  physics  involved. 

Consider  a  portion  of  an  incident  sound  wave 
scattering  from  a  finite  region  of  a  rough  surface.  The 
portion  of  the  incident  wave  might  be  a  beam,  (and  therefore 
the  beam  pattern  must  be  considered),  or  it  could  be  part  of 
a  larger  wave.  The  region  considered  is  presumed  to  be  small 
enough  so  that  the  incident  wave  may  be  regarded  as  plane, 
yet  large  enough  so  that  averaging  over  the  surface  region 
gives  meaningful  and  representative  results. 

If  the  incident  wave  ensonifies  a  large  region  of 
the  surface,  then  the  portion  being  considered  here  is  part 
of  the  integrand  of  an  integral  over  the  large  region.  This 
final  result  is  then  obtained  by  summing  or  integrating  the 
partial  results. 

Due  to  the  randomness  of  the  surface,  such  a  sum  is 
incoherent.  To  fix  our  ideas  in  the  following  discussion, 
consider  the  intensity  of  a  plane  wave  scattered  from  a 
finite  region  of  the  surface. 

There  are  two  different  contributions  to  the 
scattered  field.  One  effect  is  the  near-specular  reflection 
around  the  specular  path  (between  source  and  receiver) 
relative  to  the  mean  planar  surface.  Such  an  effect 
dominates  near  this  specular  path.  This  reflection,  often 
described  as  "glints",  is  due  to  specular  scattering  off 
large  planar  facets  of  the  surface.  When  looking  down  at  a 


sunlit  ocean  surface,  the  glints  observed  are  from  plane 
facets  momentarily  oriented  so  as  to  provide  specular 
reflection  from  the  sun  to  the  eye.  The  strength  of  this 
component  of  the  scattered  field  falls  off  rapidly  away  from 
the  specular  path. 

A  second  component  of  the  scattered  field  is 
diffuse  diffraction  due  to  the  small  scale  roughness  of  the 
surface.  This  effect  is  not  as  strongly  angle  dependent  is 
the  first  component  described  above,  and  accordingly  tends  to 
become  the  important  component  of  the  scattered  field  at 
scattering  angles  far  from  the  specular  angle. 

"Tends"  was  used  in  the  preceding  sentence  because 
as  the  surface  becomes  smooth,  and  approaches  a  plane,  the 
diffuse  component  vanishes,  while  the  "coherent  component" 
described  earlier  reduces  to  the  specular  reflection  from  a 
plane.  If,  as  suggested  earlier,  the  theory  is  describing  an 
integrand,  to  be  integrated  over  some  region,  then  this 
coherent  component  must  reduce  to  a  delta  function  within  the 
integral.  When  the  integral  is  evaluated,  the  delta  function 
produces  the  desired  specular  reflection  from  the  limiting 
plane  surface. 

Going  to  the  other  extreme,  as  the  surface  becomes 
rougher,  with  steeper  slopes,  then  the  specular  component  is 
weakened.  This  is  because  the  planar  facets  are  oriented 
over  a  larger  range  of  angles,  and  accordingly  fewer  are 
properly  oriented  for  specular  reflection  between  source  and 
receiver. 

At  the  same  time,  for  this  increasingly  rough 
surface,  the  diffuse  component  contains  a  larger  fraction  of 
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the  total  incident  energy,  and  as  the  surface  roughness 
increases,  approaches  a  scattered  field  nearly  independent  of 
angle. 

Some  little  sketches  in  Figure  1  borrowed  from 
Beckman  and  Spizzichino1 6  suggest  this  transition. 

The  first  scattering  component,  or  near  specular 
component  described  above  is  often  referred  to  in  the 
acoustical  literature  as  the  "coherent  field",  while  the 
second  is  termed  the  "incoherent  field".  Both  contribute  to 
the  scattered  average  intensity,  while  only  the  coherent 
field  contributes  to  the  mean  pressure  field,  which  pro¬ 
pagates  along  the  specular  path. 

The  acoustical  distinction  between  the  coherent  and 
incoherent  fields  is  perhaps  not  entirely  clear,  particularly 
when  discussing  the  average  intensity.  However,  there  is  a 
genuine  difference  between  the  two  fields,  because  it  is 
believed  that  different  physical  mechanisms  give  rise  to 
them. 

The  experiments  which  provide  the  most  convincing 
testimony  to  this  involve  radar  backscattering  from  the  moon. 
Using  time  delay  to  determine  which  annulus  on  the  moon  is 
illuminated  by  the  radar  pulse,  and  thereby  determining  the 
angle  of  incidence,  the  following  results  are  obtained. 

Electromagnetic  theory  predicts  that  when  a 
circularly  polarized  wave  is  reflected  from  a  plane,  the 
reflected  wave  is  also  circularly  polarized,  but  in  the 
opposite  sense.  The  radar  return  from  the  moon  however  has  a 


component  which  is  depolarized,  that  is,  circularly  polarized 
in  the  same  sense  as  the  transmitted  wave.  This  depolarized 
component  is  the  incoherent  component  or  diffuse  component, 
and  the  properly  polarized  component  is  the  coherent 
component.  A  comparison  between  the  power  in  the  polarized 
component  and  the  depolarized  component  is  shown  in  Figure 
217.  Observe  the  zero  or  minimum  delay  occurs  at  normal 
incidence,  which  for  backscattering  is  the  specular  angle. 
Note  that  the  angular  dependence  of  the  depolarized  component 
is  much  weaker  than  the  polarized,  and  exhibits  no  spike  near 
the  specular  angle,  as  does  the  polarized  component. 
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Figure  2.  A  comparison  of  the  expected  or  polarized  com¬ 
ponent  P(t)  of  the  echo  power  at_23  cm  wavelenqth 
and  the  depolarized  component  D(t).  Note  _the 
absence  of  an  initial  spike  in  the  curve  for  D(t) 
corresponding  to  the  quasi-specular  scattering 
observed  in  P(t).  (From  Hagfors,  'Radar 
Astronomy' ) 
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PROBLEM  FORMULATION  AND  BASIC  RELATIONS 


Consider  the  wave  equation  for  the  potential  U  of  a 
sound  field, 

V2u  -  —  -2-H  =  -4*Q(R,t), 

C2  at2 


where  R  »  (x,y,z),  c  is  the  sound  speed,  and  4*Q  is  the 
source  density  distribution.  The  velocity  v  and  the  pressure 
p  are  given  by 


v  *  -vu,  p 


/ 


where  p  is  the  density  of  the  medium. 
We  shall  assume 


Q(R,t)»Q(R)e-1<Jt,  U(R,t)-U(R)e 

# 

so  that  we  obtain  the  Helmholtz  equation 
v2u  +  k2U  »  -4*Q(R) , 


(4.1) 


For  simplicity  of  exposition  in  the  following,  c 
will  be  assumed  constant.  Since  the  primary  concern  will  be 
the  interaction  of  an  incident  wave  at  the  surface,  relative 
to  a  distant  source  and  receiver,  corrections  to  be  intro¬ 
duced  due  to  refraction  are  fairly  obvious.  For  example,  a 


term  representing  cylindrical  spreading  loss  to  a  point  on 
the  surface  should  be  replaced  by  the  correct  spreading  loss 
along  a  suitable  ray  path  from  the  source  to  the  surface 
position;  an  incident  angle  90  of  a  wave  arriving  at  the 
surface  should  be  related  to  the  angle  at  the  source  that  the 
corresponding  ray  makes  with  the  vertical,  the  wave  number 
k*w/c,  uses  the  sound  speed  at  the  surface,  etc. 

The  use  of  ray  terminology  already  suggests  that 
the  frequencies  involved  are  high  enough  that  a  ray 
description  of  the  sound  field  makes  sense,  which  gives  a 
lower  limit  of  100  hz,  or  perhaps  50  hz. 

The  basic  relation  governing  solutions  of  (4.1), 
derived  using  Green's  formula,  is 


U(R)  *  /  Q(R')G(R,R')dR' 
V 


+  T*  /lG<R-*>TO(r>-°(r,!K(R'c>>dr’ 


(4.2) 


where  R'eV,  the  volume  containing  the  source  Q  and  the 

a 

observation  point  R,  and  reZ,  the  boundary  surface  of  V. 
denotes  differentiation  with  respect  to  the  exterior  normal 
to  the  surface  2.  The  function  G  is  the  Green's  function,  a 
solution  of  (4.1)  when  Q»6(R-r'),  6  the  Dirac  delta  function. 

Considerable  simplification  of  (4.2)  can  be 
achieved  by  a  suitable  choice  of  the  Green's  function. 
Suppose  I  is  a  plane  surface.  Then  two  convenient  choices 


GA(R,R') 


ik  R-R' 


ik  Ri -R' 


(4.3) 


where  R1  is  the  point  obtained  by  reflecting  R  with  respect 
to  the  plane  l,  (if  l  is  the  plane  z*0,  then  Rx  =  (x,y,-z)). 

Then  it  results  that 


&G+I  , 

"StT”!  *  °'  G-|  *  °' 

I  I 

(for  l  =  {z=0},  -  -gf). 


(4.4) 


Using  either  G+  or  G_,  (4.2)  becomes 
U(R)  *  /  Q(R')G+(R,R')dR*  +  ^  /  G+(R,r)||(r)drr 

V  I 


(4.5) 


U(R)  *  /  Q(R')G-(R,R')dR' 
V 


.  3G_ 

T5t  ^  U(r)-^j-(R,r)dr, 

I 


(4.6) 


If  the  source  function  Q(R)  is  zero  outside  a 
sphere  of  radius  L,  and  if 


t*f'  f*T 


«  U 


(4.7) 


then  the  source  term  in  the  expressions  (4.2,  4.5,  4.6) 


satisfies 


ik  R-R* 


/  Q(R’)  V— 
V  R-R’ 


ik  R 


1R*  *  d ( K)- 


(4.8) 


where  d(K) 


/  Q(R' )e'lK*R,dR’ ,  K 
V 
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l|  The  region  of  space  where  the  inequalities  (4.7) 

are  satisified  is  called  the  Fraunhofer  zone  relative  to  the 
source.  Equation  (4.8)  states  that  -  at  sufficiently  large 
distances,  an  extended  source  looks  like  a  point  source. 


* 


I 


One  can  verify  directly  that  a  solution  (4.1)  with 
Q(R)«0  is  given  by  a  plane  wave,  of  the  form 

U(R)  »  Aelq'R,  (4.9) 


where  q  »  (qx,  qy,  qz),  and  q  2  «  k2. 


One  of  the  useful  properties  of  plane  waves  is  that 
a  spherical  wave,  as  in  (4.8),  can  be  represented  in  terms  of 
plane  waves.  In  fact 
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ik 


1* 


i  r. 

Tit  *  e  — a — ‘ 

— »  “7 


(4.10) 
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where  q2  ■  ^k2-qx2-qy2.  So  when  discussing  a  wave  incident 
upon  a  surface,  if  only  a  restricted  region  of  the  surface  is 
considered,  the  incident  wave  may  be  usefully  modeled  as  a 
plane  wave. 


For  reference  purposes,  let  us  describe  the 
propagation  of  a  plane  wave  in  an  infinite  space,  consisting 
of  two  half-spaces,  separated  by  the  plane  z»0,  with  density 
and  sound  speeds  p,  c,  for  z>0,  and  p2f  c2  for  z<0.  Let  U, 
02  denote  the  corresponding  potentials. 

Continuity  conditions  at  the  interface  lead  to  the 
boundary  conditions 
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3n  =  Ti T#  pU  *  p2U2 


(4.11) 


(where  ^  ^ ,  as  before). 

Let  k  =  w/R,  k2  =  t*>/c2,  K  -  ( kx  ,  ky,  kz)  = 
k(sin0cos  $ ,  sin0sin<t>,  -cos$)/  K'  *  (kx,  ky'  _,5Z)  •  K2  = 
(k2(sin02cos  4>,  sin02sin4>,  -cos02),  where  0,  02  are  the 

incident  angles  of  the  incident  and  transmitted  waves 
respectively ,  and  4  the  azimuthal  angle,  relative  to  some 
fixed  vertical  plane. 

K,  K' ,  and  K2  are  the  wave  vectors  of  the  incident 
reflected  and  transmitted  waves,  respectively  (see  Fiqure  3). 
Of  course,  all  these  vectors  lie  in  the  same  vertical  plane 
containing  the  source  and  the  receiver  at  R,  defined  by  the 
azimuthal  angle  4>. 

Snell's  law  defines  02  in  terms  of  0,  namely 

slne2  ,  Sine 

c2  c  * 

If  the  incident  wave  is  Uq  “  Ae*K  R,  then  the 
reflected  wave  is 

Ur  »  V(0)AelK'*R,  (4.12) 

and  the  transmitted  wave  is 

UT  -  W(0)AeiK2*R.  (4.13) 


-26- 


The  field  in  the  upper  half-space  is  given  by  U-  + 


Ur. 

V(0),  W{0)  are  given  by 

( p 2/p ) cos 9  -  /(c/c2)2  -  sin29 

V  =  -  -  - 

( p 2/ p ) cos 9  +  /{c/c2)  2  -  sin20 

w  =  _ 2cos9 _ 

(p2/p)cos0  +  /(c/c2)2  -  sin20 

Let  t)  *  (pc)/(p2c2)  denote  the  surface  impedance 
(the  ratio  of  the  acoustic  impedance  of  the  lower  half  space 
to  the  upper).  Then  a  perfectly  free  surface  (p2c2+0) 
corresponds  to  and  a  perfectly  rigid  surface  (p2c2-*®) 

corresponds  to  t)=0. 

The  boundary  conditions  (4.11)  for  a  perfectly  free 
surface  (p2=0)  reduce  to 

U  *  0,  (4.16) 

and  the  corresponding  reflection  coefficient  V  *  -1.  (G_ 

satisfies  this  condition.) 


(4.14) 

(4.15) 


For  a  perfectly  rigid  surface,  the  boundary  condi¬ 
tions  are 


5U  _ 
Fn 


0, 


(4.17) 


(G+ 


and  the  corresponding  reflection  coefficient  V  *  +1 
satisfies  this  condition.) 
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PERTURBATION  SOLUTION  FOR  SLIGHTLY  ROUGH  SURFACES 


A.  Perfectly  Free  Surface18 

Suppose  the  rough  surface  z  =  C(r),  r  =  (x,y)  is 
such  that  the  boundary  condition  U(r,C(r))  =  0  can  be 
approximated  by 


U(r , 0)  +  C ( r )  (r ,0)  =  0 


(5.1) 


Now  represent  U  as  the  sum  of  a  mean  field  and  a 

fluctuating  field,  U  *  U  +  u,  where  <U>  *  U,  <u>  =  0.  Also 
assume  <C>  =  0. 


Then  averaging  (5.1)  and  subtracting  the  resulting 
equation  from  (5.1)  results  in 


u(r ,0)  +  5tJ^'0)  C(r )  -  0 


(5.2) 


Subtracting  (5.2)  from  (5.1)  and  averaging  gives 


U(r,0)  +  <|H  (r,0)C(r)>  -  0 


(5.3) 


In  deriving  (5.2),  the  term  was 

neglected.  A  sufficient  condition  for  this  (i.e.,  for  juj 

»  |f?c|,  etc*>'  is 


ok  /  1 


-  (sine  +  |J)2|  <<  1 


(5.4) 
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Here  0  is  the  angle  of  incidence,  k  is  the  acoustic 
wave  number,  a  is  the  r.m.s.  surface  height,  c2  =  <C2>,  and  l 
is  the  correlation  length  of  the  surface,  as  in  (1.22). 

For  ki.  >>  2n,  this  reduces  to 

<?k  cose  <<  1  (5.5) 

For  k£  <<  2 tc ,  this  reduces  to 

ok||  ■  2nj  »  2*r  «  1  (5.6) 

where  T  -  j  -  r.m.s.  slope  of  C. 

If  sinO  21,  (near  grazing)  this  reduces  to 


ok  h  a  /IjE  «  1  (5.7) 

The  parameter  2akcos9  introduced  in  5.5  is  the 
Rayleigh  parameter. 

a)  The  Mean  Field 


Using  the  Green’s  function  G_  in  (4.4)  for  the 
plane  z*0,  the  boundary  condition  (5.3)  results  is  (r  « 
(x,y,z) ,  r*  »  (x ' ,y' )  -  (x’,y',0)). 


U(R)  "  T5  h  1  *T~T'  ~  (r')C(r')dr'  (5.8) 

1  Z  a'-O  R-r'  dz* 


Taking  the  z-derivative  of  (5.8),  multiplying  by 
C(r),  and  averaging,  the  boundary  condition  (5.2)  gives  the 
equation 


U(r) 


(5.9) 


lim  -  J  — —  ■■  —  W(p) -  (r+p)dp, 

-o  »■*  ■-«  /|7p7^ 

where  r'  =  r+p. 

Assume  the  incident  field  is  a  plane  wave,  so  that 
the  mean  field  has  the  form 

U(r,z)  =  A{elK*R  +  VeiK'*R},  (5.10) 

where  K,  K'  are  as  defined  in  Section  IV.  (The  fact  that  the 
reflected  wave  is  directed  along  K'  results  from  inserting 
(5.10)  (with  K'  replaced  by  an  abitrary  vector)  into  (5.9). 
That  is,  the  mean  field  is  directed  along  the  specular  path, 
which  is  a  result  generally  true.) 

Using  (5.10)  to  express  U  in  (5.9),  the  resulting 
equation  can  be  solved  for  the  reflection  coefficient  V, 
resulting  in  the  expression 

V  ■  -1  +  2cos0ti^(  9, $) ,  (5.11) 

where  the  "effective  admittance"  of  the  rough  surface  is 
defined  by 


e,<t>) 


a  2  +® 

a  ■*-?// 

z+0  -• 


ik/ 

P 

2+z2 

e 

4 


W(  p)elki*pdp 


2+z2 


(5.12) 


where  kj.  »  (kx,  kv) 


To  evaluate  (5.12)  for  a  general  W(p)  would  require 
detailed  knowledge  of  W.  But  if  W  depends  only  on  |p|,  i.e., 
the  surface  is  assumed  isotropic,  one  of  the  integrations  in 
(5.12)  can  be  performed,  the  z  differentiation  and  limit 
taken,  with  the  result  that 

a>  x)cy 

lye)  =  k2c2  -  ika2  J  ^  [w(y)  J0  (kysinS)  ]dy  (5.13) 

(The  derivation  of  (5.13)  depends  in  an  essential 

manner  on  W'(0)  =  0.  The  reader  should  recall  the  discussion 
on  smooth  surfaces,  and  the  corresponding  requirements  on  W. ) 

If  now  the  correlation  length  X  is  introduced  as  a 
scale  factor,  W(y)  =  w{^)#  various  limiting  forms  of  (5.13) 
are  derived  by  Bass: 

For  k A  < <  1,  or  equivalently  2n*  <<  X, 

„  ,  **,,2  _  '  t  dWiiid,  (5-)4) 

*  o  *  ax 

For  k£  >>  1,  and  (■£  -  9)  >>  — — , 

*  /Tel 

2  k2tf2cos9,  (*  ■  ^  -  0  is  the  grazing  angle)  (5.15) 

or,  recalling  5.11, 

V(0)  *  -1  +  2k2o2cos29. 


(5.16) 


-\  .  \ 


:w-v  fv^v^'*’.  >?cryy?r 


Finally,  for  k£  >>  1,  and  kl<\>2  <<  1,  (grazing 

incidence) , 


2k2 o2  ^-3«i/4  j“  J__  dW 


T)  5  ■  ■  ■  ■■  e 
^  /Juki 


0  /x 
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(x  )dx 


(5.17) 


(5.16)  will  occur  subsequently  as  a  limiting  form 
in  Section  VI  for  surfaces  with  a  small  Rayleigh  parameter. 


b) 


The  Intensity  of  the  Fluctuating  Field  for  a  Free 
Surface 


Begin  with  the  formula  (5.8)  derived  earlier  for 
the  fluctuating  component  of  the  field,  u(R),  except  now 
restrict  the  integral  to  the  actual  scattering  surface  S 
ensonified  by  the  incident  wave,  so 


.  .  1  5  f 

u(r)  *  7«  Si  I 


eik  R'r'  80 


S  | R-r • )  dz1 


(r ' ) C (r ' )dr ' 


(5.18) 


It  is  convenient  to  place  the  source  at  a  heiqht 
z0  above  the  surface,  at  (0,0,z0),  and  the  receiver  at  (D,z), 
see  Figure  4.  In  Figure  4,  R^,  R2  denote  vectors  connecting 
an  arbitrary  point  r  of  S  to  the  source  and  receiver 
respectively,  (from  the  source,  to  the  receiver).  Then 


-  j*i |  *  /z»2+M2' R! "  1*1  -  /*2+lD-rl2 


(5.19) 


-33- 


•jC*  *"«  *fc ■  ■*\  ■  « •  .  ■  »  *  %’  •,*  %*  •  * 


*:■ 


•v-;^ 

•  -  ».o 


«Y\\ 


.1 


s’.v  J 


1 


•W.y 

e-iT#1 


r  .  ,  - 


M| 


.'i 


Mv'O 

■-.’-'.'J 


.V.s  . 


■.  v 


a. 


are  the  corresponding  unit  vectors 


From  (5.11)  it  follows  that  the  reflection 
coefficient  V  differs  from  -1  by  an  amount  which  is  second 
order  in  the  perturbation  parameter  (e.g.,  the  Rayleigh 

parameter).  Accordingly  U(R)  in  (5.17)  may  be  replaced  by 
the  planar  result  for  a  perfectly  free  surface,  for  a  point 
source  at  (O,O,z0)  given  in  (4.3),  which  in  the  present 
coordinate  system  has  the  form 


U0(R)  * 


ik  R-S  ik  R-S' 

i  e 


(5.21) 


where  S  =  (0,0,z0)  is  the  location  of  the  source,  and  S'  = 
(0,0,-Zq)  is  the  image  source. 


Then  -5Y“(r,0)  *  -2ikz0  - - - ,  for  kz0  >>  1. 

Similarly,  for  kz  >>  1,  taking  the  z-derivative  in  (5.18) 

gives 


ikR  i 


k2  r  e 
u(R)  -  ~  - 


ik( R ! +R2 ) 2  z 


C(r)dr 


(5.22) 


(5.21)  immediately  implies 


<  u  ( R )  j  2  > 


//  eik(Rl+R2-Ri'“R2,) 
SxS  R12R22Hi,2R2'2 


<C(r)C(r' )>dr  dr', 

(5.23) 


where  R} ' ,  R2'  are  the  norms  of  the  vectors  connecting  the 
point  r'  in  S  to  source  and  receiver  respectively. 


Let  p  =  r'-r  be  a  new  variable.  Then  <C(r)C(r*)>  * 
c2W(p).  Now  if  the  dimensions  of  S  are  greater  than 
Xx,Xy,  (the  correlation  lengths  of  C),  the  p  integration 
may  be  extended  to  infinity,  (W(p)  ~  0  for  jp|  >>  Ax, 

Ay).  Further,  if  min(A,  v4-£+)  <<  Ri ,  R2 1  then  Ri'2  2  R^2, 

R2 ' 2  -  R22*  Finally,  expanding  the  exponent  in  terms  of  p, 

Ri  +  R2  “  Ri’  -  R2’  “  [§~  -  *  P  (5.24) 


if 


kA  2cos2e,  kA  2  «  M  (5.25) 
x  y 

where  M  »  2Ri R2/( R^+R2 ) .  Note  that  /M/kcos29  and  /M/k  are 
the  (x-y)  dimensions  of  the  Fresnel  zone  relative  to  the 
source  and  receiver.  So  (5.25)  implies  that  the  field  is 
considered  only  in  the  Fraunhofer  zone  relative  to 
irregularities  of  dimension  Ax,  Av. 


Collecting  all  the  above,  (5.22)  reduces  to 


2  a  2 


a  ^p 

2  _  k**c2  2  2 


+  a> 


<  |  u(R)  |  2>  =  — 2-  /  — - - -  dr  J  W(p)e~lk(Pl~aJL)#pdp  (5.26) 


S  Ri2R22 


(Recall  that  a 


*1  _  -2  n 

“  r7  s  (r7'  =  (“i*  az)' 


R' 


and  $=  —  = 
r2 


(lr'  "  ({V  pz)*) 


The  consideration  of  various  limiting  cases  is 
useful  in  obtaining  a  deeper  understanding  of  (5.26).  Before 
introducing  these,  recall  the  assumptions  made  in  deriving 
(5.26),  namely  (5.25),  and  the  inequalities 


L  »  *  -  max(*x,Ay);  min(*,4|I)<<  Rlf  r2«  (5.27) 

where  L  is  a  characteristic  dimension  of  the  scattering 
region  S. 

Now  suppose  kA  >>  1,  large-scale  irregularities. 

Then,  in  (5.26) 


/  dp  *  (£)* 


(•r)2 


«(|=E  -f-), 

R2  Ri 


(5.28) 


recalling  (5.20). 


Note  that  *  Pj.  defines  the  point  r  *  r0  of 
specular  reflection  in  S  relative  to  the  source  and  receiver, 
(see  Figure  5). 


In  the  opposite  case,  kJl  <<  1,  (fine  irregular¬ 

ities),  the  exponent  in  the  p-integral  is  small  for  all 
angles,  and 

+00  CO  ^>Qp 

J  W(p)e”l5c(^l_al)  *pdp  2  /  W(p)dp  =  I  t  J  W(C,h)d^dTi 

—oo  —oo  *  y  —  a> 


s  Ci  JL 
u  x  y 


(5.29) 


where  C  is  constant,  approximately  one, 


Here  the  scale  factors  Jlx,  Jty  have  been  intro¬ 
duced,  so  that  W(x,y)  »  w(j-,  ^-)  =  W(£,ti). 

x  y 

i)  "Small"  Scattering  Surface  S 

Suppose  in  addition  to  (5.25)  and  (5.27), 

k*L  cos29,  L  «  2R1R2/(R1+R2).  (5.30) 

(These  inequalities  are  somewhat  weaker  than 
requiring  that  the  entire  surface  S  be  in  the  Fraunhofer  zone 
relative  to  the  receiver.  Bass  refers  to  the  region 
described  by  (5.30)  as  the  spectral  partition  zone.) 

In  this  case,  the  r-dependence  of  the  integrand  in 
(5.26)  becomes  so  weak  that  the  r-integrand  may  be  approxi¬ 
mated  by  evaluating  the  integrand  at  an  arbitrary  point  in  S, 
and  multiplying  by|s|,  the  area  of  S. 


In  this  case,  it  is  convenient  to  reinterpret 
|u(R)  2,  and  rather  than  think  in  terms  of  a  path  connecting 
source  and  receiver,  regard  the  expression  in  (5.26)  as  a 
function  of  the  incoming  and  outgoing  wave  vectors,  in  the 
directions  given  by  o,p,  respectively.  Here  we  regard  <x,p  as 
directed  to  the  center  of  the  small  scattering  surface  S. 
Denoting  < j u ( R ) j 2  >  by  j(a,P), 
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raw 


the  result  is 


J(a,P)  =  ■ 2  |s|az2Pz2P1(k(Pi-ai)).  (5.31) 

R1  R2Z 

Fj  is  the  spatial  power  spectrum  of  the  surface ,  that  is, 
(see  1.13), 

Fi(K)  *  — — -  /  /  W(p)elK*pdp  (5.32) 

( 2n) 2  — 

The  appearance  of  the  vector  MP^-a^)  as  the 
argument  of  Fx  in  5.31  implies  that  in  first  order  perturba¬ 
tion  theory,  the  average  intensity  of  the  scattered  field 
depends  upon  only  one  component  of  the  surface  spectrum. 
This  is  a  type  of  spatial  resonance. 

In  the  special  case  kX  >>  1,  estimating  (5.32)  by 
(5.28),  (5.31)  reduces  to 

2  2 

J(a,P)  «  4*  q  Isla  2P_26(p  -a  )  (5.33) 

R12R22  1  1 

so  the  long  irregularities  in  the  surface  scatter  only  in  the 
specular  direction,  oj^  ■  Pj_. 

Conversely,  if  kX  <<  1,  from  (5.29),  the  Fj-term  in 
(5.31)  ~  a2*2.  in  this  case  the  scattering  is  diffuse,  and 
is  proportional  to  (ko)2  •  (kX)2. 

ii)  "Large"  Scattering  Surface  S 

Now  suppose  that  the  ensonified  region  S  is  so 
large  that  some  portion  of  S  lies  outside  the  spectral  parti¬ 
tion  zone,  that  is,  (5.30)  is  not  satisfied. 


In  this  case,  in  (5.26),  the  r-integral  over  S  must 
be  •  evaluated .  For  reference  purposes,  (5.26)  may  be  re¬ 
written  in  terms  of  the  power  spectrum  F1  as 


<|u(R)|2>  »  4k**  /  — — l - - — - — 


(5.34) 


A  generalization  of  this  result,  for  composite 
surfaces,  will  be  encountered  subsequently. 


Suppose  kX  >>  1,  the  case  of  large  irregularities. 
Then,  using  the  estimate  (5.28), 


u(R)  |  2>  =  4k2a2  /  2  —  -  ~)dr 

1  S  RX2R22  r2  1 


Changing  coordinates  in  the  integral,  (see  Appendix 
A  for  details),  gives  the  result 


<|u(R)|2>  -  =  (l*°g>.se>2. 

•  I  o  2/«  .2  ,  ~  ^ 


Pq»2 (^1 0+^20 ) 2  (Rl0+^20)2 


(5.35) 


R0  in 
2o  ^ 
R1  0 


Here  a0,  p0, 
D-r0 

S,  where  -s -  - 

«2  0 

cose . 


Ri0,  R2o  refer  to  the  specular  point 
ro  .  2o  . 

—0  *  «'  r°  *  Also  V  - 


Observe  that  in  (5.35)  the  term  in  the  denominator 
which  gives  the  spreading  loss,  between  source  and  receiver  is 
no  longer  Ri2R22  (compare  (5.31)),  but  has  switched  to 
(Ri  +  R2)*,  which  gives  the  smaller  spreading  loss  observed 
in  specular  reflection  from  a  plane. 
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Also  note  the  occurrence  of  the  Rayleigh  parameter 
2kacos0  (assumed  small,  recall)  in  the  numerator  of  (5.35). 
Observe  that  here  Lambert's  law  (scattering  ~  cos';)  is 
applicable. 


In  the  opposite  limit,  kl  <<  1,  fine  irregular¬ 

ities,  (5.28)  implies  that 


<  j  u(R) | 2> 


C(k*  )(k*  ) (ka)2  a  23  2 

- * - 1 -  /  -5 — —  dr 

n2  S  R12R22 


(5.36) 


Clearly  in  this  case,  the  scattering  is  very 
diffuse,  with  the  incident  angle  0  of  little  importance,  and 
the  entire  region  S  influencing  the  average  intensity. 

B)  Rigid  Surface19 

Now  suppose  the  initial  condition  is 


^  (r,C(r) )  »  0, 


(5.37) 


where  the  normal  derivative  d/3n  is  given  by 


(5.38) 


,+V*Y 


Here  y_  * 


Again  represent  U  *  U  +  u,  the  mean  field  plus  a 
fluctuating  field.  Assuming  S(r)  is  suitably  small, 
expanding  (5.37)  in  C,  averaging  and  differencing,  the 

boundary  conditions  (to  first  order)  for  0,  u  become 


au(r,0)  _ 

5z  ' 


v  u>  - 
r  &z2 


(5.39) 


du 


(r  ,0) 


=  Y*7rU  - 


(5.40) 


a) 


Here  T  -  V  *  (V  V'  ’r  *  ly)- 

The  Mean  Field  for  a  Rigid  Surface 


Using  the  Green's  function  G+  in  (4.4)  for  the  mean 
plane  z-0,  the  boundary  condition  (5.40)  results  in 


u(R) 


— r-7'  {c(r')^Hi£li  -  y(r ' )  V  ,  U(  r ' )  }dr ' . 
R-r' 

(5.41) 


This  expression  for  u(R)  can  be  used  in  conjunction 
with  the  boundary  condition  (5.39)  to  obtain  an  integral 

equation  for  U  for  z«0,  analogous  to  (5.9)  for  the  case  of  a 
free  surface.  It  is  not  particularly  revealing,  however,  and 
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will  be  omitted.  Nevertheless,  if  one  assumes  an  incident 
plane  wave,  with  a  reflected  wave,  as  in  (5.10),  it  can  be 
shown  that  Snell's  law  continues  to  hold,  that  is,  the 

reflected  component  of  the  mean  field  U  for  an  incident  plane 
wave  propagates  along  the  specular  path. 

If  the  surface  is  assumed  isotropic,  i.e.,  w(p) 
depends  only  on  jp  j,  the  situation  simplifies  to  the  point 
where  a  formula  for  the  reflection  coefficient  V(0)  can  be  >; 

obtained,  namely,  (compare  with  5.11), 

■ 

i  % 

a 

(5.42) 

.  <r 
.  f 
.V 


(9  is  the  specular  angle  between  source  and  receiver.) 

The  effective  impedance  ng  is  given  by 
2  ®  ik 

Ti  (9)  *  k2a2  -  I  {-k3sin9cos9j1  (kysin9)W(y) 

9  *  o  y 

+  ^^-[(ik3y-k2sin20)Jo(y)  +  ^lll£(  i-iky-k2y2 ) ^  (kysin0)  ]}dy. 

(5.43) 

As  before,  letting  the  correlation  length  X  be  a 
scale  factor,  W(y)  *  w(^)*  various  possiblities  exist: 

For  kX  <<  1 


0  >> 


For  kl  >>  1 ,  4»  =  - 


/TcT 


^g(e)  5  k2o2cos39, 


or,  referring  to  (5.42) 


V(  9 ) 


1-k2g2cos26 

1+k2a2cos29 


(5.45) 


(5.46) 


Finally,  for  kA  >>  1,  and  near  grazing,  <|//JcX  <<  1, 


ng(6) 


a2/ET  e  J  _1 dW( x )dx 

2JtVZw  0  x3/2  dx 


(5.47) 


A  difference  between  the  representation  of  the  mean 
field  for  a  free  surface  and  a  rigid  surface  occurs  here 
because  the  reflection  coefficient  V(9)  for  a  rigid  surface 
can  have  a  pole  near  grazing,  namely  (see  5.42), 


cose  »  -hg(0)  (5.48) 

Since  jTlg(®)j  <<  1  in  the  cases  considered  here, 

a  good  approximation  to  a  solution  of  (5.48)  is  given  by 
setting  tig  *  0,  and  iterating,  with  the  result  that  the 
angle  9p  at  which  the  pole  occurs  is  given  by 


The  effect  of  the  pole  on  U  is  as  follows:  For  a 
point  source  above  the  boundary  and  at  large  distances,  the 


reflected  component  of  U  has  the  form,  for  9  *  0  , 

r 

ikR  n 

U  (R)  =  V(  0  )  +  0  (-— ) .  (5.50) 

r  R  kR 

* 

where  R®»  /x  j  2  +y !  2  + ( z+z  0 ) 2  ,  and  9  is  the  specular  angle 
between  the  source  (at  (O,O,z0))  and  receiver  (at  (x,y,z)), 
see  Figure  5. 

But  if  V  has  a  pole  (when  920p)  ,  then  the 
reflected  component  has  the  form 


Ur(R) 


(i+y}) 


iit-s2 


where 


s 


2 


(?)  [  1  + 


“g® 


(5.51) 


(5.52) 


The  following  asymptotic  formulas  are  valid  for 
grazing  incidence,  <1**0: 


*  Note  that  S  -  (1  +  |-)Rio 


(1  + 


20  • 


See  (5.35)  and 


s 


<<  1 


(5.53) 


So  if  R,  the  position  of  the  receiver,  is  such  that 
the  specular  angle  020p,  and  if  js  j  is  large,  the  reflected 
mean  field  will  experience  additional  attenuation. 

However,  8p  <  -j,  and  as  grazing  incidence 
approaches,  0  <  6  s  the  reflection  coefficient  approaches 
1,  the  attenuation  due  to  the  pole  disappears,  and  (5.50)  is 
the  appropriate  estimate. 

In  fact,  as  will  be  shown  in  the  next  section,  the 
average  intensity  for  a  rigid  surface  equals  that  for  a  free 
surface,  at  grazing  incidence. 

b)  The  Intensity  of  the  Fluctuating  Field  for  a  Rigid 

Surface 


The  relevant  formula  for  u(R)  is  (5.41),  except  now 
restrict  the  integral  to  the  actual  scattering  surface  S,  so 


.  _ik  )R-r ' 

«<»>  -Tt  i  Vrr 

s  JR-r1  j 


{ C  ( r  * ) 


. , a2u(r ' ) 


dz’2 


-  Y(r,)Vr,0(r')}dr'. 


(5.54) 


Introduce  the  same  geometry  as  before,  (see  Figure 
4),  with  Rx ,  R2 ,  «,  p  defined  in  (5.19),  (5.20).  As  in  the 
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case  of  a  free  surface,  U  may  be  replaced  by  the  zero-order 
aproximation  CJ0,  where 


ik  |R-S  j  ik  |R-S*  j 

uo(R)  3  ”Tr-3_i|  +  {r-s '  1 


(5.55) 


as  before  (compare  5.21),  where  S  *  (0,0, z0)  *  -S’,  provided 
|cos9+rig  j  >  0,  so  that  Vh-1  is  a  valid  approximation. 


Then  (5.54)  can  be  rewritten  in  the  form 


u(R) 


0-VPj 

R1R2 


eik(R1+R2) 


C(r)dr. 


(5.56) 


But  (5.56)  immediately  implies  that 


<|u(R) | 2> 


k4q2  r 
x2  S 


Rl2*22 


dr 


7  W(p)e"lk( Pl“ai) *pdp. 

•  flD 

(5.57) 


The  assumptions  made  to  derive  (5.57)  are  the  some  as  in  the 
free  surface  case,  (5.25)  and  (5.27). 


But  (5.56)  has  the  same  form  as  (5.26),  with 
(azPz)2  replaced  by  (1-ai*Pi)2,  so  estimates  cor¬ 
responding  to  (5.31)  through  (5.36)  are  valid  in  the  case  of 
a  rigid  surface  also,  with  the  above  substitution. 


Note  that  (5.35)  is  valid  for  rigid  surfaces  also, 
because  at  r“r0,  i.e.,  o^-p^  ®ozPoz  *  1-a0i»p0i,  as  shown 


at  the  end  of  Appendix  A.  So  as  remarked  earlier,  if  the 
scattering  area  is  large  and  kJl  >>  1,  Lambert’s  law  (5.35)  is 
applicable  to  both  a  rigid  and  a  free  surface  (for  a  small 
Rayleigh  parameter,  of  course). 

C)  Surface  with  an  Impedance  Boundary  Condition19 


Assume  the  initial  condition  is 


au 

a7j(r,C(r) )  =  -ikriU(r,C(r) ) 


(5.58) 


where  n  is  the  impedance  at  the  interface  between  the  upper 

& 

and  lower  half  spaces,  and  ■g-jj  is  the  normal  derivative 
defined  in  (5.38). 


Letting  U  ■  U+u  as  before,  relatively  few  changes 

occur.  The  boundary  condition  for  U,  5.39,  has  a  term  -iknU 
added  on  the  right  side.  The  effective  impedance  rig  is 
changed  by  the  addition  of  the  impedance  t{  on  the  right  side 
of  (5.43). 

And  finally,  the  expression  a228z2  for  a  free 
surface,  in  (5.26),  or  (1-aj.iPx)2  in  (5.57)  is  replaced 
by  a  more  complicated  expression  j  F(  a,P,k,k2,*,Jl2  j  2  given  in 
Appendix  B.  So  the  expressions  in  (5.31)  through  (5.36)  are 
valid,  with  a22Pz2  replaced  by  j F J  2 .  Note  that  in  this  case, 
(5.35)  becomes 


<ju(R) 


4k*o* J  P( O0  ,  P0)  ■ 
Pqz2(rio+r2o) 2 


(5.59) 
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VI 


SCATTERING  FROM  SURFACES  WITH  LARGE-SCALE  ROUGHNESS 


The  Kirchoff  method  or  tangent  plane  solution20. 

To  discuss  scattering  from  a  rough  surface  when 
there  is  no  small  parameter  available,  the  following  approxi¬ 
mation  is  useful. 

Assume  the  surface  is  so  smooth  that  at  each  point 
r  of  the  surface  S  the  wave  field  U(r)  may  be  represented  as 
the  sum  of  the  incident  field  U0(r)  and  the  field  reflected 
from  the  tangent  plane  at  r,  that  is, 

5U(r) 

U(r )  »  (1+V)U0(r),  2gLLL  *  (l-VJ-^fr)  (6.1) 

Here  is  the  normal  derivative,  and  V  is  the 
reflection  coefficient. 

If  one  assumes  further  that  the  incident  wave  is  a 
plane  wave  at  each  point  r  of  S,  with  a  wave  vector  K 
depending  on  r,  (6.1)  becomes 

U(r)  »  (1+V(6))U0(r),  .  i(K*N)( l-V(6))U0(r)  (6.2) 

where  N  denotes  the  normal  to  S  at  r,  6  is  the  angle  between 
K  and  N,  and  the  plane  wave  approximation  is  used  for  U0,  so 
that  the  normal  derivative  s  iK *N  U0. 

The  above  representation  requires  that 
1/3 


cos 9  >>  (ka)~ 


(6.3) 


where  a  is  the  local  radius  of  curvature  of  S  at  r.  This 
condition  results  from  the  requirement  that  the  tangent  plane 
at  r  is  a  good  approximation  to  S  over  a  region  whose 
dimensions  are  large  relative  to  the  wavelength  of  the 
incident  wave. 


Let  the  scattering  surface  S  be  complemented  by  a 
plane  S'  and  the  hemisphere  Cr',  (see  Figure  6)  so  that  a 
closed  region  is  constructed  containing  all  the  field 
sources.  Then  if  R  is  an  arbitrary  point  within  the  reqion, 
by  Green's  theorem  the  field  U(R)  has  the  form 


U(R) 


-  U0(R) 


I 

S+S'+CR 


(6.4) 


But  since  all  the  sources  are  within  the  region,  it 
follows  that 


Tn  S+S^Col^0^^(Vi)  -  £nUo(r)ldr  *  °* 


(6.5) 


Subtracting  €.5  from  €.4,  the  result  is 


(6.6) 


The  Kirchoff  approximation  consists  in  assuming 
that  on  S1,  U«Un,  so  the  integral  over  S'  is  zero. 

Further,  assuming  a  radiation  condition,  that  is, 

gikR ' 

at  infinity  both  U  and  U0~ - ,  R'eCp,  the  integral  over 

R' 

Cri  may  be  neglected.  Using  6.2  to  represent  U-U0  on  S, 
gives  the  result 


U(R) 


ikR 

uo<R)  +  in  J  U0(r))dr, 

s  R2 


(6.7) 


where  V(r)  is  used  to  indicate  the  dependence  of  the  reflec¬ 
tion  coefficient  on  the  local  angle  of  incidence.  Here  r  is 
the  length  of  r'  along  S,.as  seen  in  Fig.  7. 

If  the  source  is  taken  as  a  point  source,  then  6.7 
may  be  written 


U(R)  *  U0(R>  +  jz  J  VU)^  [--k[Rr--2)]dr 

S  Rj  R2 


(6.8) 


where  (see  Figure  7),  R^ 2  =  ( 2o-C(r ) ) 2+J r j 2,  R22  *  (z-C(r))2 
+  jD-r  j2. 

Now  assuming  kRj,  kR2  >>  1,  (the  term  used  to 
describe  this  is  that  the  source  and  receiver  are  in  the  wave 
zone  relative  to  S),  the  integrand  in  (6.8)  can  be  approxi¬ 
mated  by  -i(N*q)exp(ik(Ri+R2)  )/(RiR2)  where  N  is  the  normal 
to  S  at  r,  as  shown  for  example  in  Fig.  6,  and  is  given  by 


A)  The  Mean  Field  U 

a)  "Small"  Scattering  Area 

Suppose  the  dimensions  L  of  the  surface  S  are  less 
than  the  dimensions  of  the  Fresnel  zone,  that  is 

kL2cos26  <<  r1#  kL2cos20 '  <<  R2.  (6.13) 


Then  in  (6.12),  k(Ri+R2)  2  k(R10+R20)  -  qi  •  r,  where  Rio* 
R20  are  the  distances  between  a  fixed  position  on  S  and  the 
source  and  receiver,  respectively,  and  now  q  is  a  constant 
vector,  over  all  of  S.  (6.12)  then  becomes 


_ieiJc(Rl0+R2°)  .  ,  . 

U(a'P) - Tur77r2-q -  l  v(r)(qz-qi-Y)exp(-'i(qi-r+q2C(r))dr. 

s0 

(6.14) 


Here  U(a,P)  rather  than  U(R)  is  used  to  in  indicate  that 
(6.14)  represents  the  field  in  the  direction  p  after  S  is 
illuminated  by  an  incident  plane  wave  with  wave  vector  ka. 


But  as  <y>  «  0  and  C  and  y  are  independent,  it 
follows  directly  that 


U(a,p) 


-iq  veiIc*Rlo+R2o) 


l 71 R  l  qR2o 


where  fi  is  the  first  order  characteristic  function  defined 
in  (1.6).  But  the  integral  above  equals  (2u)25(qi),  and 
qi  ■  0  only  for  the  specular  path  connecting  source  and 
receiver,  so  finally 
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-iuq  V(eo)elk(Rlo+R2o) 


4  0*20 


■6(q±)e 


-2k2c2cos20( 


(6.15) 


where  60  is  the  angle  of  specular  reflection  at  q^  =  0.  In 
(6.15),  fi(<3z)  was  evaluated  for  a  Gaussian  surface,  as 
qz  =  2kcos0o  when  qi  a  0. 


(6.15)  again  asserts  that  the  mean  field  propagates 
along  the  specular  path.  Further,  the  effective  reflection 
coefficient  is  V( 0o)exp(-2k2a2cos20o) ,  so  that  the  coherent 
component  of  the  scattered  field  decays  exponentially  with  an 
increase  in  the  Rayleigh  parameter.  Note  that  (5.16)  is  a 
special  case  of  (6.15)  for  small  Rayleigh  parameter,  for  0 
away  from  grazing,  and  kA  >>  1. 


'Large"  Scattering  Area 


Now  suppose  S  is  so  large  that  the  inequalities  in 
(6.13)  are  not  satisfied.  In  this  case,  to  discuss  the  mean 
field,  return  to  equation  (6.9)/  rather  than  continue  from 
(6.12). 


Assuming  z0,  z>C,  the  spherical  waves  in  (6.8)  may 
be  expanded  in  plane  waves,  namely. 


eikRx  •  A  +-  ei(KT+tz0-C(r)]K2) 

~T~m  7*  L - iC - dK 


(6.16) 


ikR2  i  +-  i(K»(D-r)+fz-C(r) JK  ) 

—  - «: - J-dK 


(6.17) 


where  K  -  (Kx,Ky),  and  Kz2  -  k2-  {K  |2. 
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in  6.9 


As  before,  the  normal  derivative  has  the  form  given 


Using  6.9,  6.16  and  6.17  in  6.8,  and  projecting  the 
integration  from  S  onto  the  mean  plane  S0(z=0),  the  result  is 

i  r  r  +“  [Kz+K'z-Y*(K-K')] 

U(R)  -  U0(R)  +  -i—  /  [V(r )  JJ— - jy-rv-r- - 

16*3  S0  —  Vz 

exp{ i [ ( K-K ' )  *r  +  K*‘D  +  z0K„  +  zK„ ' -C(r ) ( K  +K  • )  ]  }dKdK' ]dr . 

2  Z  2  2 

(6.18) 

To  calculate  ti  from  (6.18),  again  use  that  <y>  *  0, 

and  C  and  y  are  independent  random  variables.  For 

simplicity,  set  V®  1 .  The  average  of  exp(-i (K  +K„ ' )  ( C (r ) )  = 

z  z 

fi(K,+K,’)  as  before,  and  the  integral  of  exp ( i ( K-K ' ) • r ) 
z  z 

gives  (2x) 26(K-K' ) .  The  final  result  is 


i  fi<2Kz)  r  , 

U(R)  =  UQ  ( R)  r  .Oxp[i(K*D+(z0+z)K2idK 


(6.19) 


Recalling  that  fj(2Kz)  =  exp(-2k2a2cos29)  defines 
the  effective  reflection  coefficient  in  the  previous  section, 
(6.19)  asserts  that  the  mean  field  is  the  sum  of  the  incident 
field  U0(R)  and  the  superposition  of  reflected  plane  waves, 
arriving  from  different  angles,  and  with  their  own  reflection 
coefficient. 


Let  R'  «  / I D| 2+( z0+z ) 2  be  the  distance  from  the 
mirror  source  (0,0,-zq)  to  the  receiver  at  (D,z).  Then  it 
has  been  shown  that  (6.19)  reduces  to 
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J»(R)  -  O0(R)  +  - 


(6.20) 


XJ(R)  =  O0(R)  +  - - exp(-2k^cos^90) , 

R' 

where  90  is  the  specular  angle  connecting  source  and 
receiver,  if 


2  (ka)4sin26o+(ka)^(3cos0o-1) I  <<  kR'. 


(6.21) 


For  (ka)2  small,  (6.21)  reduces  to 


«  1/ 


(6.22) 


which  is  a  very  weak  restriction  on  a. 

However,  for  (ka)2  >>  1,and  for  9o  bounded  away 

from  0  and  x/2,  that  is. 


en  » 


ka/7 


and  \  -  0O  « 


(6.23) 


(6.21)  is  equivalent  to 


2(ka)lfsin220o  «  kR’, 


(6.24) 


which  is  a  good  deal  more  stringent  than  (6.22). 

This  is  due  to  the  fact  that  for  (ka)2  large,  the 
effective  reflection  coefficient  inside  the  integral  in 
(6.19)  will  vary  rapidly  over  S  unless  (6.24)  is  satisfied. 
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B)  The  Intensity  of  the  Fluctuating  Field 

Returning  to  the  representation  of  U(R)  given  in 
(6.12),  note  that  for  large  Rayleigh  parameter,  here  written 
qza ,  the  main  contribution  to  the  integral  will  come  from 
the  points  of  stationary  phase,  namely  the  points  for  which 
Vr[k(R1+R2)-qzC(r)]  =  0,  that  is 

qzY  +  qx  -  0,  (6.25) 

recalling  that  y  =  VrC,  and  a j_  =  7rRi»  01  *  -VrR2. 

But  the  solution  to  (6.25),  To  *  ~qi/qz  defines 
a  point  of  specular  reflection  between  the  source  and 
receiver.  So  in  (6.12),  the  term  V(r)  (qz-qi*Y)  may  be 
evaluated  for  y  a  solution  of  (6.25).  Since  qz-qj.*Yo  * 
|q|2/qz»  (6.12)  becomes 

U(R)  *  J  V(R|[]qr) l2exp(irk(R1^R2)-q2C(r)] )dr,  (6.26) 

where  V(r)  now  denotes  the  reflection  coefficient  evaluated 
at  the  angle  of  specular  reflection  at  the  point  r. 

But  from  (6.26),  the  following  formula  for  the 
averaged  intensity  u  results. 


<  u ( R)  2>  «  <  U(R)  2>  -  <U( R) >  2 


«  — ~J  /  ■  — *'-"k  I2  fo'  I2  exp[ik(R1^R2-Ri R2'3*  (6.27) 

(4,1)2  S(J  Sq  RlRl'R2R2'qsq8» 

[  <exp(-iq  C(r  )+iq  •  C ( r  * )  )>  -  <e”iqzC (r }  ><eiqz '  C  {r ' )  >  ]drdr ' . 


t 


For  large  Rayleigh  parameter,  (qzo)^  >>  1 
expanding  W(p)  in  a  Taylor  series,  the  result  is  that  (6.29) 
is  valid  for 

kRr  2  >>  i,  kRcos2er  2  >>  i  (6.3D 

x  y 

where  R  2  R1 ,  R2,  cos29  2  cos20'.* 

Interpreting  the  inequalities  in  (6.31)  geometri¬ 
cally,  the  result  is  that  for  (qzc)2  >>  1,  (6.29)  is  valid 
if  the  region  where  specular  reflections  dominate  substan¬ 
tially  exceeds  in  size  the  first  Fresnel  zone. 

If  now  we  consider  the  case  (qza)2  «  1,  the 
expression  in  (6.30)  reduces  to  qz2cr2w{p),  and  in  this  case 
(6.29)  nearly  coincides  with  (5.26).  A  factor  |q  |4/16  occurs 
in  (6.29),  and  a  term  k**az2Pz2  in  (5.26).  These  agree 
when  ® 1 ,  «Z*-PZ  (scattering  in  the  specular  direc¬ 
tion),  but  not  elsewhere.  Further  the  conditions  for 
validity  of  (6.29)  turn  out  to  be  (5.25)  and  (5.27),  that  is, 
the  region  where  specular  reflections  dominate  must  be  much 
smaller  than  the  Fresnel  zone. 

Not  only  is  (6.29)  valid  at  reasonable  ranges  for 
both  large  and  small  values  of  the  Rayleigh  parameter,  but  as 
will  be  shown  below,  (6.29)  reduces  to  the  correct  limiting 
value  as  the  mean  slope  of  S  approaches  zero.  Further,  for 
large  values  of  the  Rayleigh  parameter,  a  limiting  form  of 
(6.29)  can  be  derived  involving  the  slope  probability  density 
function  which  agrees  with  results  derived  using  an  alterna¬ 
tive  approach  to  scattering. 


For  all  of  these  reasons,  (6.29)  is  a  fundamental 
result  of  scattering  theory. 

a)  Small  Scattering  Surface 

Now  suppose  that  the  dimensions  of  the  scattering 
surface  S,  and  its  projection  Sq  are  sufficiently  small  that 
the  r-dependence  in  the  integral  in  (6.29)  may  be  neglected. 
Then  (6.29)  becomes 

(V  I2  |S0  |  |q  |4  +»  _iq  ,p  2i 

*<■•»>  -  7«*~«"K2q  ~  --  'I'**.’-*.”'-'*  ‘V*’ 

2  (6.32) 

Here  as  in  Section  V,  the  notation  J(a,0)  is  used 
as  a  reminder  that  the  average  intensity  is  now  a  function  of 
the  incoming  and  outgoing  wave  vectors,  a  (in),  p  (out),  and 
q  *  k(p-a).  |Sq  j  is  the  area  of  S0. 

For  (qzo)2  »  1,  (6.32)  is  valid  if 


L  COS0  L  COS0'  L 

- ,  -*  R” -  «  r^.  RlCOS0,  R2COS0’  «  ^ 


where  Lx,  Ly  are  x-y  dimensions  of  S0 


For  (q zo)2  «  1,  (6.32)  is  valid  if 


(6.33) 


kL  l  cos 2  0  kL  X 


(6.34) 


plus  similar  inequalities  with  Ri,cos@  replaced  by  R2,cos0'. 
(6.34)  corresponds  to  the  restriction  (5.30),  derived  earlier 
for  the  perturbation  solution. 
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Now  consider  (6.32)  for  (q2c)2  >>  1.  Then  the 
terra  f x 2 ( q2 )  »  exp(-qz2c2)  can  be  neglected.  Further,  as 
soon  as  1-W(p)  becomes  significantly  different  from  zero,  the 
f2  term  (see  (6.30))  is  also  essentially  zero.  So  for  large 
qza,  only  the  region  near  p=0  needs  to  be  considered. 


But  then,  taking  a  Taylor  expansion  of  1-W(p) 
around  p»0,  and  evaluating  the  integral  for  large  (q2a), 
using  the  method  of  stationary  phase,  (or  actually,  Laplace's 
method  here),  the  result  is 


J(«,P) 


V!2|S°1  hi  “ _ 2. 

(4*RiR2q2)2  <322°2/|Wxx{0)Wyy(0) 


exp[- 


2q. 


2c2 


(6.35) 


But  setting  W  (0)  »  W  (0)  ■  ,  and  r  2 

ax  ^  2  ll  |  2  A 

2  -  ‘  *  7 


a2/l  2 ,  r  2  a  a2/«2  (g.35)  can  be  rewritten  as 

xy  y 


J(a,P) 


so|M 


4R12R22qs(‘* 


-w2(Yq)/ 


(6.36) 


where  w2(y)  is  the  slope  probability  density  function  defined 
in  (1.2),  and  Yo  *  ~qi/qz  *s  the  facet  slope  providing 
specular  reflection  between  source  and  receiver. 


The  derivation  of  (6.36)  requires  that 


!? 

>; 


0 

K\ 


t_3 


jM 


) 


in  (6.36)  is  independent 
(6.36)  expresses  the  high 
frequency  limit  of  the  integral  in  (6.32). 

Radar  backscattering  studies  from  the  moon  exhibit 
some  frequency  dependence,  even  near  the  specular  direc¬ 

tion21,  so  to  that  extent  (6.36)  is  not  correct. 

Further,  the  derivation  of  (6.36)  given  requires 

very  little  information  regarding  W(p),  only  that  second 
derivatives  exist.  (But  see  an  alternate  derivation  of 
(6.36)  below.)  So  if  the  scattering  depends  upon  more  than 
three  parameters  a,  Ax,  Ay,  then  (6.36)  is  too  simple. 

However,  its  very  simplicity  is  appealing. 
Further,  the  fact  that  the  slope  density  function  w2  is 
evaluated  at  y0  *  -<3i/<Jz  is  also  interesting.  A  facet 

with  slope  y q  is  such  that  the  wave  vectors  a  (from  the 

source),  and  p  (to  the  receiver)  make  equal  angles  with  the 
normal  to  the  facet.  So  the  reflection  from  such  a  plane  is 
specular  reflection  from  source  to  receiver.  Thus  (6.36)  can 
be  interpreted  as  the  product  of  a  geometric  term  Jq  |Vqz4, 
an  area  term,  a  reflection  coefficient  and  a  spreading  loss 
term,  all  multiplied  by  the  probability  that  there  is  a  facet 


2a)  >  1  if  q 


hi 


<  1 


and 


o  ,  hi  L  ki  I 
<V”2  »  lf  57- 


Here  r  *  max(Tx,  Ty). 


Note  that  the  expression 
of  k.  This  is  due  to  the  fact  that 
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suitably  oriented  to  provide  specular  reflection  between 
source  and  receiver. 


This  is  an  expression  of  the  idea  that  for  a  large 
Rayleigh  parameter,  the  scattering  is  due  to  glints  from 
properly  oriented  facets  of  the  surface.  And  since  the  mean 
plane  is  flat,  <C>  *  0,  this  probability  is  highest  near  the 
mean  specular  path,  so  this  term  should  dominate  there. 

(6.36)  is  of  sufficient  interest  that  an  alterna¬ 
tive  derivation  will  now  be  given. 

Begin  with  U(a,P)  as  given  in  (6.14). 

Then 


|U(a,P)|2 


1 

MitR^)2 


I  I  |v|2(q2-q1-Y)(qz-q1-y1)- 
s0  So  (6 

exp[-iqi(r-r' )+q2(C(r)-C(r' ) ) ]drdr' 


Now  in  the  exponent,  let  C(r)  -  C(r')  2  VrC(r-r')  *  -y*p, 
where  p=*r'-r  as  before.  Neglecting  the  p  dependence  in  y* 
results  in  a  p  integral  of  the  form 


(2x) 26(qi+qzY) r 


(6.39) 


again  assuming  that  the  integral  is  dominated  by  the  contri¬ 
bution  near  p*0,  so  the  integral  may  be  extended  to  the 
infinite  plane. 

Introducing  (6.39)  into  (6.38),  and  then  averaging 

with  respect  to  the  random  variable  y,  the  result  is 

/ 


J(<x,0)  =»  <|[J(a,|3)  |  2>  =  - - - -/|v|  2(qz-qi*Y)26(qi+qzY)w2(Y)dY 

4R^ ‘R2 

V( e ) | 2  j S0  j  jq 
4R12R22qz4 

where  Yo  *  ~qi/qz  as  before.  90  is  the  specular  angle 

relative  to  the  normal  to  the  facet  with  slope  Yo*  that  is, 
N  •  Q 

cos  9  o  ■  "7q|*  In  use  was  made  of  the  fact  that 

(qz-qi*Yo)2  *  |q|4/qz2*  The  estimates  required  in  deriving 

(6.40)  are  precisely  the  inequalities  (6.37)  obtained  in 
deriving  (6.36). 


-w2(yq)  / 


(6.40) 


Note  that  in  deriving  (6.40),  no  mention  was  made 
of  the  correlation  function  W.  However,  the  underlying 
premise  throughout  the  whole  section  is  that  the  surface  S  is 
smooth  enough  that  the  tangent  plane  approximation  makes 
sense.  But  as  remarked  earlier,  a  reasonable  smooth  surface 
requires  W  to  have  second  derivatives.  So  the  derivation 
leading  to  (6.40)  is  not  more  general  than  that  for  (6.36). 

b)  Shadowing  for  Small  Scattering  Surface22 

In  deriving  the  above  relations,  such  as  the 
formula  for  jO(o,P)  j2  in  (6.38),  the  implicit  assumption  was 
made  that  the  entire  surface  S  was  illuminated,  so  that  the 
integration  could  be  taken  over  the  entire  projected  area 
Sq  • 

However,  for  rough  surfaces,  with  steep  slopes,  and 
near  grazing,  obviously  some  portions  of  the  surface  are  in 
shadow,  and  do  not  interact  with  the  incident  wave.  To 


correct  this  deficiency,  a  factor  relating  to  the  probability 
of  being  shadowed  should  be  introduced. 


The  simplest  way  to  introduce  the  shadowing  cor¬ 
rection  is  to  assume  no  diffraction  and  modify  (6.38)  so  that 
the  integration  is  now  over  S^xS.^,  rather  than  S0xS0, 
where  is  the  projection  of  the  illuminated  portion  of  S. 
Continuing  as  above,  approximating  C(r)-C(r')  s  -yp,  and 
using  (6.39),  we  arrive  at  (6.40).  But  now,  rather  than 
computing  the  average  with  respect  to  y  using  the  density 
function  w2(y),  the  y  average  is  computed  using  the  effective 
slope  density  function  we(y;a,P)  which  is  the  probability 
of  the  occurrence  of  a  facet  with  slope  y  for  which  an 
incoming  ray  a  is  not  intercepted  by  some  other  portion  of 
the  surface,  and.  for  which  the  outgoing  ray  p  is  likewise 
unshadowed. 

Now  the  result  is 

we(y;a,P)  -  w2(y)P(y;a,p)  (6.41) 

where  P(y;a,P)  is  the  probability  that  neither  of  the  rays  a 
or  P  reaching  the  facet  with  slope  y  intersect  the  surface 
elsewhere. 


The  only  case  where  theory  has  produced  an  evalua¬ 
tion  of  P(y;a,P)  which  has  been  verified  experimentally  is 
when  o,  p  lie  in  the  same  vertical  plane,  the  plane  of 
incidence23. 

Let  o  ■  (sin0cos4,  sin0sin4,  -cos0),  and  p  ■ 
(sin0'cos4',  sin0'sin4',  cos0')r  where  0,  0'  are  the 

incidence  angle  of  a,  p  with  respect  to  the  normal  to  the 
mean  plane,  and  4  ,4'  are  the  angles  a,  p  make  with  respect 
to  the  x-axis,  (see  Figure  3). 


If  <t>*  *  it+$  (backscattering )  and  ‘5’  =  max(0,0')» 

P(Y?®#P)  *  Htcot'S-Y  cos4>-Y„sin4>)  • - 3 - .  (6.41) 

X  y  1+A("S,  4>) 


If  4*'  * 
P(y;«#P) 


4>  (forward  scattering) 

=  H(cot9,-Yvcos4>-Yvsin<t>*  )H(cot0+Y  cos4>+Y  sin<t>) 
x  y  x  y 

•  - - -  (6.42) 

1+A  (  0 , 4>  )+A(  0  '  ,  $  1  ) 


Here  H(x)  is  the  unit  step  function,  1  for  x>0,  0 
for  x<0.  A(0,<|>),  for  a  Gaussian  surface,  is  defined  by 


pf  *3  ^ 

A(0,4>)  “  e  2  -aErfc(~) ] . 


(6.43) 


where  a2  ®  cot29/<rx2cos24>+ry2sin24>)  • 

The  difference  between  the  two  cases  is  as  follows. 
In  (6.41),  (backscattering),  if  the  lower  ray  is  not 
shadowed,  then  certainly  neither  is  the  upper  ray,  that  is, 
the  two  events  {a  not  shadowed}  and  {p  not  shadowed}  are 
related,  one  is  a  subset  of  the  other. 


However,  in  forward  scattering,  (6.42),  the 
assumption  is  made  that  the  two  events  are  independent,  which 
seems  reasonable  in  this  case. 


No  shadowing  theory  has  been  developed  for  an 
arbitrary  relation  between  ♦  and  4'.  The  difficulty  lies  in 
making  the  transition  between  two  independent  shadowing 
events,  and  the  totally  dependent  arrangement  in 
backscattering . 
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In  the  absence  of  any  analysis  of  the  transition# 
the  suggestion  has  been  made  to  use  (6.42),  treating  the 
events  as  independent,  as  long  as  p  lies  outside  some  wedge 
containing  a,  i.e.,  as  long  as  j  >  e0  >  024.  Presumably 
e 0  need  not  be  very  large,  e.g.,  n/12;  although  Bass  suggests 
e0  *  -j.  For  the  reverse  inequality,  (6.41)  should  be  used. 

One  further  note,  in  the  final  form  of  (6.40), 
where  w2(yo)  is  replaced  by  we(Yo»a>P)*  the  step  functions 
in  (6.41)  and  (6.42)  are  one,  as  the  arguments  are  positive. 

But  when  these  results  are  extended  to  a  large 
scattering  surface,  these  step  functions  will  play  a  role. 

The  shadowing  theory  presented  thus  far  is  based 
upon  the  approximation  in  (6.38)  of  C(r)-C(r')  by  -y*P.  This 
eliminates  the  random  variable  C  in  favor  of  y-  And  the 
introduction  of  a  correction  for  shadowing  is  reasonably  well 
understood,  when  averaging  over  the  slope  vector  y. 

However,  as  remarked  earlier,  the  formula  (6.32), 
(or  (6.29)),  seems  more  fundamental  than  (6.40),  in  that 
(6.32)  will  give  either  (6.40)  or  the  perturbation  result 

(5.31) ,  according  to  the  size  of  qza.  But  to  the  author's 
knowledge,  no  shadowing  correction  has  been  applied  to 

(6.32) . 

There  is  a  formula  for  we(C;a,p)  analogous  to 
(6.41),  where  we(C?<*,0)  is  the  probability  density  function 
of  the  surface  height,  given  that  the  rays  o,p  are  not 
shadowed,  namely. 


where 


;a,P)  *  (l+A(0,0)+A(0*, ♦•))•[}  w(  C  '  )dC  '  ]  ^A(  9  '  *)+A(  9  '  '  *  '  ]  ^ 


+“ 

•  //P(Y;a,P)dr 


(6.45) 


The  fact  that  P(C;<*,P)  depends  upon  C  will  prevent 
the  averaging  of  the  term  C(r)-C(r')  giving  f2,  as  before  in 
(6.32). 

One  suggestion  for  introducing  a  shadowing  correc¬ 
tion  into  (6.32)  or  (6.29)  is  to  approximate  an  integral  over 
sill  as  follows: 

<  J  {  }dr>  =  /  P(«,P)  <{  }>dr, 

Sill  S0 

where  P(a,p)  is  the  probability  of  the  rays  a,  p  reaching  an 
illuminated  portion  of  the  surface  independently  of  the 
height  C,  and  slope  y. 


For  an  isotropic  Gaussian  surface,  and  for  a,p 
independently  shadowed. 


P(a,P) 


COt  9  ' 

/  w2(Y)dY 
-cotQ _ 

1 +A  { 9 ) +A{ 0 ’ ) 


(6.46) 


If  the  shadowing  of  a  and  the  shadowing  of  P  are 
dependent,  then 


P(a,P) 


9 


(6.47) 


cot¥ 

/  w2 ( Y )d y 


l+Af'S) 


where  "S'  =  max (9,9’)  as  before. 


The  shadow  corrected  version  of  (6.32)  is  then 


J(a,p)  *  J0(afp)P(a,p), 


(6.48) 


where  J0(a,p)  is  given  in  (6.32). 


For  a  non-isotropic  surface,  the  integrals  in  the 
numerator  in  (6.46)  and  (6.47)  are  replaced  by  the  integral 
of  W2(7x'Yy)  restricted  by  the  step  functions  appearing 
in  (6.42)  and  (6.41),  respectively. 


Large  Scattering  Surface 


If  the  inequalities  (6.34)  are  not  valid,  so  that  a 
large  region  is  involved  in  the  scattering,  then  the  r 
integral  needs  to  be  evaluated.  The  result  is 


< | u(R) | 2>  »  / 

So 


J(«/P) 


(6.49) 


where  J(a,P)  is  given  in  (6.32),  or  for  large  (q2a) ,  in 
(6.40),  where  no  shadowing  correction  is  made;  or  by  (6.40) 
multiplied  by  P(y;<*,P),  or  (6.48),  which  is  (6.32)  multiplied 
by  P(«,P). 


It  is  interesting  to  consider  the  limit  of 
< | u(R)  2>  as  r2+0  (but  letting  kR-*-®  so  (6.31)  is  still 
satisfied).  Let  us  write  (6.49)  in  the  form 
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t  s*  v' 


£5 


JW 


L-  I  P(«>P)|v(r)[i|aJJl  +J  e-^l*p[f2(qz,-qz;P) 


< | u ( R )  1  2 >  L- 

'  '  (4n)2  S0  (RiR2qz)2 

-  fi2(qz) ]dPdr 


(6.50) 


As  r2->-0,  f2  (see  (6.30))  becomes  sharply  peaked  at 


p=0,  so 


J  e”iql*P[f2(qz,-qz;p)  -  fi2(qz)]dp  s  ( 2*) 26 (q± ) [ 1-f x 2 (qz ) ] 

■oOO 

(6.51) 

But  as  in  (5.35)  qi®0  occurs  at  the  point  of 
specular  reflection  with  respect  to  the  mean  plane ,  namely  at 
,  z°  . 

c°  * 

The  result  is 


<  u(R) 


4|v(60)|2(1-f12(q2)) 


(6.52) 


( R* o+^20 ) 


where  we  have  used  the  fact  that  at  qi^O,  |q|=qz*p2. 

zo 

0O,  as  before,  is  the  specular  angle  at  r0,  cos0o  *  -g —  * 
z  *  ® 

R2  o  • 

When  (qza)2  >>  fi250,  and  (6.52)  gives  the 

correct  limiting  result  in  specular  reflection  from  a  plane. 


In  (6.52),  the  probability  P(a,P)  +  1  as  r2-*-®,  as 
A(9)-*-0  as  r2-*-®. 
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VII. 


EXTENSIONS  TO  THE  THEORY:  COMPOSITE  SURFACES^3 


There  are  basically  two  solutions  to  the  rough 
surface  scattering  problem  that  have  gained  acceptance.  One 
is  the  perturbation  solution  (5.31),  for  a  small  Rayleigh 
parameter,  and  the  other  is  (6.40),  for  a  large  Rayleigh 
parameter.,  using  the  tangent  plane  approximation.  (Actually 
(6.40)  should  be  replaced  by  the  shadowed  version,  with  w2(y) 
replaced  by  we( y ; a, p ) . ) 

As  has  been  remarked  several  times  already,  both 

(5.31)  and  (6.40)  can  be  obtained  from  (6.32).  (At  least  the 
free  surface  perturbation  solution  can  be  obtained  from 

(6.32) ,  to  within  a  geometric  multiplier,  although  since  the 
rigid  surface  solution  coincides  with  the  free  surface 
solution  at  grazing  incidence,  perhaps  the  difference  is 
slight. ) 

However,  (6.32)  has  not  been  used  extensively  as  a 
scattering  solution,  probably  because  a  Fourier  transform 
must  be  evaluated,  whereas  (5.31)  and  (6.40)  are  simple 
function  evaluations. 

Now  (6.40  -  6.44),  and  (5.31)  are  valid  in  dif¬ 
ferent,  nearly  complementary  regions  of  the  parameter  space 
(k,a,p,o),  namely  the  geometrical  optics  solution  is  valid 
for  large  ko,  near  the  specular  angle,  while  the  perturbation 
solution  is  valid  for  small  ka,  and/or  near  grazing,  away 
from  the  specular  angle. 

Clearly  it  would  be  desirable  to  have  a  simple 
solution  which  would  bridge  the  gap  between  the  two  existing 


solutions,  as  the  result  would  be  a  theory  valid  over  all 
parameter  space. 

Note  that  the  two  solutions  are  responding  to 
different  features  of  the  scattering  surface.  The  geometri¬ 
cal  optics  solution  is  describing  specular  reflection  from 
properly  oriented  planar  facets  of  the  surface,  while  the 
perturbation  solution  expresses  the  diffraction  from  smaller 
features  of  the  surface. 

This  notion  led  to  the  introduction  of  a  composite 
surface,  that  is,  expressing  the  surface  as  a  sum  of  a  large 
and  a  small  surface. 


2  =  C(r)  =  CL(r)  +  Cs(r) 


(7.1) 


where  in  (7.1)  and  subsequently,  the  subscripts  L  and  s  refer 
to  the  large  and  small  surfaces  respectively. 

The  idea  here  is  the  large  surface  contains  the 
planar  facets,  (long  waves  in  surface),  and  superimposed  on 
these  are  the  small  ripples  (short  waves)  responsible  for 
diffraction. 

Again  note  the  utility  of  assuming  a  Gaussian 
surface,  for  if  both  Cl  and  Cs  are  Gaussian,  so  is  C. 


The  first  person  to  use  this  notion  of  composite 
surface  was  Kuryanov26.  He  suggested  modifying  the  small 
surface  (perturbation)  solution  in  two  ways:  i)  Express  the 
geometric  factor  in  the  solution  relative  to  the  normal  to 
the  large  surface  (Nl  *  N(yl))*  rather  than  in  terms  of 


I-.-' 


~  % 


the  vertical  normal  to  the  mean  plane;  and  ii)  since  the 
orientation  of  the  large  facets  is  random,  average  the 
modified  perturbation  solution  in  (i)  with  respect  to  y. 

Kuryanov  wrote  before  an  adequate  shadowing  theory 
was  developed,  which  limited  his  results.  However,  Brown2 
recently  independently  rederived  Kuryanov's  results, 
including  shadowing. 

An  interesting  way  to  develop  Brown's  results  is 
due  to  Dashen3.  Begin  with  (6.32),  repeated  here  for 
reference. 


V  z  Sn  q  H  +«  •  ?2 

J(a,P)  *  -  /e"1^!  P[exp(-q  o2(1-W(p))-e"qz  a  ]dp. 

(4uRlR2q  )2 

Z  (7.2) 


Now  (7.1)  implies 

a2W(p)  *  °L2wL( P)+°s2ws( P) • 


Assuming  (qzas)2  is  small,  then 
exp(q22ag2Ws(p))  5  1+<3z2as2ws (  P )  • 


(7.3) 


(7.4) 


2  2 

Then,  neglecting  the  term  e”qz  °  in  (7.2),  we  can 


write 


,5,  .  ill  s°  q  k  [*j  .-Hi-o.Wd- 


(4nR1R2q_) 


(7.5) 


+  q  2  1  •-19i-P(e^«V<1-,fL«P»).a  ip)dp] 
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But  the  second  term  in  (7.3),  the  Fourier  transform 
of  the  product,  becomes  the  convolution  of  the  transform  of 
the  two  factors. 

Now  suppose  that  (qzc)2  is  large  enough  so  that, 
as  in  (6.35), 


7  W  -  i!i  w  ,(r/q  ;o,P) 

«  2  eL  z 
Mz 


(7.4) 


where  a  shadowing  correction  has  been  introduced  by  using 
weL'  defined  in  (6.41),  rather  than  w2l. 


Using  (7.4)  in  the  convolution  form  of  (7.3),  the 


result  is 


|v|2|s0||q|* 

J(®,P)  *  - - - - - - — w  r(Y0fa,p) 

Rl^V  eL 


v  |Sp||q 

4Ri2R22 


(7.5) 


/  FL  (q  { Y-Y0) )w_T ( Y?a,P)dy, 


recalling  that  Yo  =  -c3i/qz'  and  in  the  last  integral  the 
dummy  variable  in  the  convolution  has  been  set  equal  to 


On  comparing  the  last  term  in  (7.5)  with  (5.31),  it 
is  clear  that  this  term  is  a  shadow  corrected  averaged 
version  of  (5.31).  In  fact,  if  kA  >>  1,  then,  (just  as  in 


the  transition  from  (5.31)  to  (5.33)),  weL(yf«»P)  5  6(0), 
(as  rL2  +  0  also),  and  Fls(-qzY0)  2  6(qx).  Thus  this 
term  in  (7.5)  reduces  to  (5.33)  also. 

Clearly  this  last  term  in  (7.5)  incorporates 
Kuryanov's  notion  of  averaging  the  small  surface  result  with 
respect  to  the  slopes  of  the  large  surface.  However,  the 
geometric  term  jq  |4  is  not  included  in  the  averaging. 

Using  a  different  approach.  Brown  included  the 
geometric  term  in  the  convolution.  If  (7.5)  was  modified  by 
taking  the  |q  |4  term  within  the  y-integral,  written  in  the 
form 

i«r-<w2*o‘-OMW.  (7.6, 

this  gives  one  interpretation  of  expressing  the  geometric 
factor  in  terms  of  the  slope  y  of  the  facet  of  the  large 
surface. 

A  more  fundamental  version  of  (7.5)  could  be 
provided  by  replacing  the  term  |q  |4  by  the  correct  expression 
for  the  perturbation  solution,  i.e.,  az2pz2  for  a  free 
surface,  (1-a±*p^)2  for  a  rigid  surface,  or  more 
generally  |p(a,P)|2.  However,  in  line  with  Kuryanov's  notion 
of  expressing  these  factors  in  terms  of  the  normal  to  the 
facet  with  slope  y,  let  az,  pz  be  interpreted  as  the 
component  along  the  normal  N(y),  and  a^,  Pi  the  perpen¬ 
dicular  component. 

That  is. 


N  V* 


Recall  that  N(y)  *  (1»-y)/*|y[  +1. 

As  remarked  earlier  after  (6.31),  |q  |4  should  be 
replaced  by  16k4az28z2.  Accordingly,  the  second  term  in 

(7.5)  has  the  form 

4k4  }V  j2  jS0  j  +» 

- - - - —  /  az(7)2Pz(Y)2Flg(qz(Y-Ya) )w2l(y)  dY-  (7.8) 

Rl  R  2 

This  is  essentially  the  form  used  by  Andreeva  et 
al29,  and  before  that  Bachmann30,  and  Kuryanov26.  Observe 
that  in  the  integral  we£,  has  been  replaced  by  w2l,  as 
earlier  Andreeva  et  al31  showed  that  at  low  grazing  angles 
the  shadowing  correction  is  not  important,  except  possibly  at 
very  high  frequencies.  Since  the  modification  and  averaging 
of  jq  |4  over  y  is  only  a  significant  correction  at  low 
grazing  angles,  it  is  consistent  to  replace  w6l  by  w2I< 
above.  Neither  Bachmann  or  Kuryanov  had  a  significant  amount 
of  data  with  which  to  test  their  model.  In  reference  29, 
however,  scattering  data  for  grazing  angles  from  3*  to  10* 
was  available,  at  several  frequencies  and  wind  speeds,  and 
the  corrections  introduced  in  (7.8)  (for  the  case  of  back- 
scattering)  produced  a  model  which  fit  the  data  rather  well. 

Probably  the  various  versions  of  the  composite 
surface  solution,  ((7.5),  or  the  substitution  suggested  in 

(7.6) ,  or  (7.7)),  are  all  essentially  equivalent,  as 
weL(Y*'a'P)  for  moderate  values  of  rL2  will  be  con¬ 
centrated  near  y*0,  where  the  above  expressions  differ  very 
little. 
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One  must  observe,  as  does  Brown,  that  the  integral 
term  in  (7.5)  should  include  the  step  functions  (see  (6.41) 
or  (6.42))  in  the  definition  of  w6l,  in  order  to  be 
rigorous.  Other  than  complicating  the  evaluation  of  the 
convolution,  it  is  not  clear  whether  the  correction  so  intro¬ 
duced  is  signficant.  Some  numerical  experience  would  be 
useful  here. 


A  more  interesting  question  relates  to  (7.1),  the 
division  of  the  surface  into  large  and  small. 

The  basic  assumption  here  of  course  is  that 
(k a)2  »  1,  for  otherwise  the  perturbation  solution  would 

probably  be  valid. 

Given  this,  the  division  of  large  and  small  is 
given  by  dividing  the  spectrum  at  a  value  k*  to  be  chosen. 

FlL(K)  *  Fi<K)H(k*-|K|)  (7.9) 

and 

Flg(K)  »  Fi(K)H( | k| -k* )  (7.10) 

Recall  H  is  the  unit  step  function,  one  for  posi¬ 
tive  argument,  zero  for  negative. 


Then 


+<a 

°L2  *  !  plL(K)dK 


(7.11) 


+® 

-  /  FU(K)dK 


a 


2 


(7.12) 


Given  a  smooth  enough  spectrum  around  qj_,  further 
terms  in  the  expansion  begun  in  (7.14)  can  easily  be 
supplied,  via  Laplace's  method. 

So  if  |qj_  j  >  k*,  and  F!  smooth  near  qj_,  (7.5) 
can  be  approximated  by  function  evaluations,  eliminating  the 
need  to  evaluate  integrals.  As  |qj.  |-*-k*,  some  attention 
has  to  be  paid  to  the  approaching  cutoff  in  Fls. 
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SCATTERING  PROM  A  MOVING  ROUGH  SURFACE 
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Now  consider  a  surface  varying  in  time,  2  *  C(r,t). 
The  basic  assumption  is  that  the  variations  in  the  surface 
C(r,t)  and  slopes  y(r,t)  are  small  relative  to  the  incident 
frequency  co0  and  incident  wave  velocity  c  ■  u>0/k  respec¬ 
tively.  Here  the  incident  wave  U0  ~  e”lu)ot. 

Such  an  assumption  allows  an  immediate  generaliza¬ 
tion  of  (6.14),  namely 


,  ik(R10+R2o)-i“»ot  , 

-  4.r10-r- -  /(V( 

exp(-i(qi*r+q2C(r,t) )dr 


(8.1) 


As  remarked  several  times  before,  such  a  formula 
will  produce  correct  results  for  large  or  small  Rayleigh 
parameters,  and  will  be  the  only  form  considered  in  this 
section. 


Again,  the  dominant  contribution  to  the  integral 
occurs  near  the  stationary  point,  Yo  ■  and  the 

term  qz-qi*y  reduces  to  |q|2/qz  as  before.  Here  the 
time  variations  of  y(r,t)  are  assumed  to  be  much  slower  than 
that  of  the  exponential  exp( -iqzC( r ,t ) ) ,  so  that  the  time 
dependence  in  y  may  be  neglected.  For  a  small  Rayleigh 
parameter,  the  term  qz-qi*Y  -  |q|2/qz  and  may  be  taken 

outside  the  integral  if  kA  >>  1.  This  is  because  in  this 
case,  that  of  a  gently  sloping  surface,  the  field  is  con¬ 
centrated  near  the  specular  direction. 
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The  function  of  interest  in  the  case  of  a  moving 
surface  is  the  power  spectrum  of  the  field,  defined  by 


1  h 

S(a,P;u)  -  /  <  0(a,P;t+x)U  (<x,p;t)  >  e1<i,XdT. 

(8.2) 

Introducing  (8.1)  into  (8.2),  using  the  approxima¬ 
tions  indicated,  and  proceeding  as  in  Section  VI  in  deriving 
(6.32),  the  result  is 

| v| 2  I  So  I M “  +7  .  .  . 

S(o,P;«*>)  =  - -  /  exp  ( -xq  .  •  P  +  x  (  u)-a)0  )  x  ) 

-2ti(  4nR1R2q,)  2  “» 

z 

*  f2(q2*-qz;p»T)<3pdT  (8.3) 

where  f 2 ( qz , -q2 ; P , t )  =  exp ( -q2 2  a  2 (  1 -W ( p , t ) ) )  for  a 

Gaussian  surface. 


The  result  in  (8.3)  could  be  shadowed  by  multi¬ 
plying  by  P(a,P)  as  in  earlier  sections.  This  is  consistent 
with  the  time  dependence  of  C  and  y,  for  the  stationary 
assumption  implies  that  the  probability  density  function  for 
C,  and  therefore  for  y,  is  independent  of  time. 


Suppose  (qzcr)  <<  Then  expanding  f2  in  powers 

of  (qzo)2W,  the  result  is 


i 

S(oc,P;w) 


| 


co  q  4n 

I  -2—  S  (afp,«) 
n=0  nl 


(8.4) 


where  the  term  multiplying  the  sum  is  evident  from  (8.3),  and 


(8.5) 


S  (a,P;u>)  *  J  exp(-iqi  *p+i  {  w-io0  )  x)  (  a2W{  p,  t)  )  dpdx 

—  CD 


But  as  the  Fourier  transform  of  a2w(p,x)  is  F(K,w), 
as  given  in  (1.11)  or  (1.15),  (8.5)  asserts  that  Sn  is  the 
n-fold  convolution  of  F.  The  first  few  terms  are  as  follows: 


S0(a,P;a>)  »  (2it)  36(qi)  6(  w-cl'0)  (8.6) 

S1(a,P;u)  =  Fi  (qx)  6 ( w-a)0-to+ ( qx ) )  (8.7) 

S2(a,p;u>)  =  !— -  /  F1(K)F1(qi-K)6(aj+(K)+u)+(qi-K)-(a)-u)0)  )dK 

{2%)  (8.8) 

where  K  *  (kx,  ky),  as  before. 


Sg  represents  the  spike  directed  along  the  specular 
path,  at  frequency  u>0.  For  the  moving  ocean  surface,  if 
Fj(K)  is  assumed  zero  in  some  circle  (interval)  around  K*0, 
depending  upon,  the  wind  speed  above  the  ocean,  then  Sj  will 
also  be  zero  in  some  interval  around  <*>0.  s2  is  in  general 
not  zero  for  wswg,  and  will  fill  in  the  gap  around  the 
central  spike  left  by  Si28. 


Now  suppose  (qzff)2  >>  1.  Then,  as  before,  f2 
decays  so  rapidly  for  W(p)*1,  (p*0)  that  the  integral  in 

(8.3)  may  be  approximated  by  replacing  1  -  W(p,t)  by  a  Taylor 
series  in  p  and  t.  The  result  is  that  S  is  Gaussian  in 
nature. 


— -  exp( 


( ( W-Wg ) -«m) 2 


0 


S(a,P;u) 


1 


(8.9) 


where 


Di 

7"' 


(Aw)2  =  -2qz2o2§  . 


(8.10) 


Here  D  is  the  determinant  of  the  Hessian  matrix  of  W  (with 

respect  to  p  ,  p  ,  and  t)  evaluated  at  p=0,  t=0,  and  d  is  the 
x  y 

2x2  principal  minor  of  D.  Di  is  D,  with  the  first  row  of  D 
replaced  by  (qj./  0).  The  second  derivatives  of  W  needed 

above  can  be  calculated  in  terms  of  Fi  and  w+,  see  (1.12)  or 
(1.16). 

If  the  scattering  surface  is  not  small,  so  that  the 
r  integral  must  be  evaluated,  the  result  is 

S(R,w)  «  let  /  S(a,8;o>)dr  (8.11) 

lS°l  S0 

Again,  to  obtain  the  correct  normalization, 

particularly  if  S0  is  so  large  so  as  to  include  scattering 
near  grazing  incidence,  (8.11)  should  probably  include  a 
shadowing  correction. 

Furthermore,  if  C  is  decomposed  into  the  sum  of  a 

large  and  small  surface,  a  composite  surface  form  of  (8.3) 

could  be  derived,  which  would  combine  the  contributions  to  S 
from  the  large  surface  (using  (8.9)  evaluated  with  respect  to 
C  £,),  and  the  contribution  to  3  from  the  small  surface, 
(using  (8.4),  evaluated  in  terms  of  Cs),  perhaps  convolved 
with  (8.9).  In  particular  for  low  grazing  angles,  the  term 
Fi(<3l)  in  (8.7)  would  be  replaced  by  an  integral  as  in 

(7.8). 
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F/G  2@/l 


NL 


‘2  0 


Prom  (A. 11)  and  noting,  from  A. 5  and  A. 6 
B,  A. 8  becomes 


1 


1 


(z+z0)2  Poz2(Rio+R2o) 


which  is  the  desired  result. 


Using  the  relations  in  (A. 9)  and  (A. 10), 

2  (z0+z)2 


1-ct  ,  •  6_  .  =  1-  Ml 


*ol  ^ol 


B‘ 


B‘ 


a  *8 
02  KOZ 


that  R^o  + 

(A. 12) 


(A. 13) 


so  that  the  rigid  surface  kernel  and  the  free  surface  kernel 
agree  at  r0. 
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The  function  F(a,P)  has  the  following  form: 


2k2  F <a,p) 


k{p2-P)^1-|P1 


p/k22-k2  2  +  kp2/l-"  p  i2 


i  k2p2-k22p  r  v 

i+v>[ — pT=p —  *2vfJ  -  tv-vKV-^KI2?- 


Here  k2,  p2  refer  to  medium  below  the  interface. 
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developed  Phillips  spectrum/ 

Fx ( k )  *  B/k3,  g/u2  <  k  <  ks  (one  dimension). 

Then 

o  *  u2  / -  ■  .0049u2  ra. 

2g2 

k  u2 

2r2  2  B  log  (-|— ) 

Here  u  ■  wind  velocity,  (m/sec),  g  «  9.81  m/sec2,  B  ■  .0046 
is  a  dimensionless  constant,  ks  is  the  cut-off  between 
gravity  and  capillary  waves,  k8  2  2w/* 3  m~ 1 


Table  1 


u( knots) 

u( m/sec) 

o(«) 

r 

Jt(m) 

10 

5.1 

.13 

.096 

1.35 

20 

10.3 

.52 

.112 

4.65 

30 

15.4 

1.2 

.120 

10.0 

40 

20.6 

2.1 

.125 

16.8 

Let  k  -  2*f/c,  f  in  herts,  c  s  1500  m/sec. 
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